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_ SUMMARY 

unrf 

The basic theory of space-charge flow is described, and three gen- 
eral methods of solution are discussed. Specific space- charge -flow ge- 
ometries are described in detail; these include flow between coaxial 
cylinders, concentric spheres, coaxial cones, hyperbolic flow, flow be- 
tween inclined planes, and circular flow. Since all the solutions dis- 
cussed are for space- charge-limited flow, a discussion of partial space- 
charge flow is presented. The Pierce method of electrode design is de- 
scribed, and some practical limitations of this method are discussed. 
Electrode shapes for various space-charge flows are shown. 


INTRODUCTION 

The performance and design requirements of electrostatic rocket 
engines have brought a new emphasis to the theory of space- charge flow. 
Considerable reductions in engine size and substantial improvement in 
engine efficiency may be gained by the use of ion and charged-particle 
accelerators operating with current densities at or near the space- 
charge limit. Such high current densities bring about definite space- 
charge effects, so that Poisson T s equation must be used to describe the 
potential distribution in the engine accelerator. In addition to the 
requirement of high current density, the electrostatic rocket engine 
must be capable of operating for very long times, of the order of years 
for interplanetary space missions. Such long- duration operation re- 
quires near-perfect "ion optics" to keep electrode sputtering erosion 
to reasonable limits. Since the ions, or charged particles, are accel- 
erated to kinetic energies in the range from a few thousand to hundreds 
of thousands of electron volts, each ion striking an electrode may 
sputter away several electrode atoms. For these reasons, some estimates 
of the degree of perfection of ion optics place an upper limit of ion 
impingement of less than 0.01 to 0.1 percent of the total ion beam to 
keep erosion to a reasonable level. 
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The space -charge -flow problem in electrostatic rocket engines is 
characterized by the necessity for the ions, or charged particles, to 
leave the engine permanently through an open exhaust aperture, as shown 
in sketch (a): 



Since the ion source is at a potential considerably higher (or lower, 
for negative ions) than the potential of free space, the equipotential 
lines in the accelerator will tend to bow, or bulge, out of the exhaust 
aperture. This aperture effect not only reduces the current carrying 
capacity of the accelerator, but also introduces another complicating 
boundary condition on the space-charge-flow equations. 

The complexities of the space-charge flow in electrostatic rocket 
engines indicate that complete analytical design of space -charge -flow 
accelerators will be very difficult. For this reason, recourse has been 
made to analog or empirical methods (i.e., experimental n cut- and- try fT ) . 

Not only are these methods inadequate for extrapolation or scaling of 
engine designs, but in addition they do not offer sufficient basic under- 
standing of space-charge flow to allow positive prediction and evaluation 
of the effects of configuration changes in design optimization studies. 

It is evident that a convenient, universal method for the analytic 
solution of the space- charge-flow problem would be of great value in the 
design of electrostatic rocket engines. The survey of existing literature 
on space- charge flow reported herein has yielded no completely universal 
analytic method, but some progress has been made in that direction. 
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The first section of this report is devoted to a discussion of the 
three methods of solution that have "been used to date for space- charge- 
flow problems. The remainder of the section on space -charge -flow theory 
is a compilation of analytic solutions for specific space- charge -f low 
geometries found in the literature and also includes a calculation of 
the flow between coaxial cones that was not previously available. 

Since the analytic solutions to space- charge flow usually assume 
that the flow fills all space, it is necessary to take a segment of this 
flow and to provide electrodes external to the ion beam to force the flow 
to conform to the theoretical pattern. The second section of this report 
is devoted to the problem of electrode design. The method of Pierce is 
described and is applied to a number of space -charge -flow geometries to 
obtain electrode shapes. Some of these shapes are accurately recomputed 
using techniques found in the literature, and others appear not to have 
been computed previously. Some practical limitations to the application 
of space- charge flow theory and the Pierce method of electrode design 
are discussed. 

The parameters of the various space-charge flows examined in this 
report may be used to determine the theoretical performance of electro- 
static rocket engines. An analysis is made in reference 1 in which the 
various space- charge flows are compared in terms of relative engine per- 
formance. 

The work described in this report was carried out as part of the 
electrostatic propulsion research program of the NASA Lewis Research 
Center. 


SPACE- CHARGE-FLOW THEORY 

In the study of the performance of ion engines the analysis of 
space -charge flow is very important. Since the early studies of Child 
(1911, ref. 2) and Langmuir (1913, ref. 3), who solved the problem of 
space-charge-limited flow between parallel plates more than half a cen- 
tury ago, there have appeared extensive publications on space- charge ef- 
fects in certain new geometries. Most of the early analyses (refs. 2 
to 5) have been concerned with one -dimensional cases, and it has been 
only in recent years that two- and three-dimensional cases were solved 
or solutions were suggested (refs. 6 to 13). 

In most ion engines, ions move as a very nearly continuous fluid in 
an electrostatic field in which the magnetic force is negligible. The 
flow of ions is assumed to form a system of curves (streamlines) filling 
a portion of space such that the velocity is a single-valued function of 
position. 


4 


The fundamental differential equations governing steady-state space- 
charge flow are: 

Poisson T s equation: 

v 2 ® = 

O 

Equation of motion: 

4 

mv = q V® 

Conservation of charge: 

V j =0 

(Note that 0 = cp Q - cp. ) 

It is mathematically convenient to assume that the flow is space- 
charge limited. In this case the velocity and electric field will "both 
he zero on the emitting surface. Under these conditions it can readily 
he shown (ref. 9) that the flow will he irrotationalj that is, 


(1) 

( 2 ) 
(3) 


V X v = 0 (4) 

In this case, the velocity can be expressed as the gradient of a 
scalar action function: 


v = (5) 

Combining equation (5) with the vector identity v = (v • v)v yields 


v = (vW • v)\^ 

Application of equation (4) to the identity 

(tW • v)W = V(\^) 2 - X (V X W) 


( 6 ) 


gives 


v = -| V(W) 2 


(7) 
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Substitution for v from equation (2) yields 

V(W) 2 = (Ip) TO (8) 

Hence, 

(vw) 2 = (9) 

where the constant of integration is included by letting $ = cp Q - cp. 

The five equations used to formulate space-charge-flow problems are: 
Hamilton-Jacobi equation: 


Poisson ! s equation: 


Total ion energy: 


(w) 2 = 





(9) 


( 1 ) 


( 10 ) 


Conservation of charge: 

V * 0 = 0 

Vector current density: 

J = pv 


(3) 

(ID 


Analytical Methods 

Historically the first analytical solutions of the space- charge- flow 
problem were for particular cases of rectilinear flow between specified 
equipotential surfaces. These are flow between infinite planes (refs. 2 
and 3), flow between coaxial cylinders (refs. 3 and 4), and flow between 
concentric spheres (ref. 5). One of the first attempts at a general so- 
lution of the space- charge-flow equations was made by Spangenberg (ref. 6) 
following his formulation of the general equations in terms of the action 
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function. His multidimensional analysis yields a partial differential 
equation in terms of the action function, which is much too unwieldy 
for exact solution. 

In 1949 Meltzer (ref. 7) introduced a unique approach to the solu- 
tion of space-charge-flow problems. By his method the space-charge- 
flow parameters are expressed in terms of the acceleration and velocity. 
Meltzer shows that there are two conditions necessary for a space-charge 
solution to exist; namely, 


V X v = 0 

(12) 

<1 

<.j, 

<jl* 

II 

O 

(13) 


Using the relation v = (v • v)v to obtain the acceleration from a pre- 
viously assumed velocity, v and v are tested in equations (12) 
and (13). If they satisfy (12) and (13), they define a type of space- 
charge flow whose parameters are as follows: 


- ©■»<" • T) 

(14) 

- ©j 

(15) 

(!)# * • d? 

(16) 

■ (f) e ^ (v ■ 7) 

(17) 


This method leaves one big equation unanswered. How may a velocity 
function be found that will satisfy equations (12) and (13)? The search 
for an answer to this question led Meltzer (ref. 11) and later Rosenblatt 
(ref. 12) to develop the trajectory action function method, which is more 
readily applicable to space-charge-flow problems than any previous 
method. 

In the action function method of Meltzer and Rosenblatt it is as- 
sumed that the flow of charged particles is laminar, and that an orthog- 
onal curvilinear coordinate system 


£(x,y,z) 


(18a) 


T)(x,y,z) 
£(x,y, z) 
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(18b) 

(18c) 


can be found such that the ion trajectories are completely described by 
variation of the parameters t) and £ (i.e., q = const, and t, = const, 
specify a particular trajectory). Then the orthogonal surfaces 
i = constant will be surfaces of constant action, and the Hamilton- 
Jacobi equation (9) becomes 


0 - 


(19) 


where h^ is the scale factor corresponding to the coordinate £. Rosen- 
blatt shows that if a function T is defined such that 

T = h^4> ( 2 °) 


the general differential equation for space-charge flow can be written 

li/2 [i|( 0 f) + Br ] mZM) (2: 


where 


C 


h 2 h 3 


hn 


( 22 ) 



(23) 


K(ti,0 



(24) 


The differential equations for several known space -charge -flow cases are 
derivable from equations (21), (22), (23), and (24). The variables f, 
■n, and the scale factors, and equations (21) to (24) are tabulated 
for certain of these cases in table I. It is still necessary, however, 
to apply the classical techniques to obtain final solutions. 
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Particular Solutions of Space-Charge Flow 

It is of interest here to discuss particular solutions of the differ- 
ential equations listed in table I and to express these solutions in 
nomographic form. In addition to these previously listed space-charge- 
flow solutions, there are two other known analytic space-charge flows in- 
cluded in the following sections. 

Rectilinear flow between infinite plane electrodes . - Sketch (b) 
illustrates rectilinear flow between plane electrodes. 



In 1911, Child (ref. 2) published the solution to this case. From his 
analysis he obtained what is now known as the Child-Langmuir space - 
charge law: 


j 



(25) 


where cp Q is the potential at x = 0. When the constants are written: 


4 [ 2 ^ 

9 € °? m 


5.467X10" 8 A" 1 / 2 
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where A is the molecular weight per ionic charge, equation (25) "be- 
comes 

i = 5.467X10" 8 A _1 / 2 ( C p o _ cp) 3 / 2 x " 2 (26) 

A nomogram relating j, (cp Q - cp). A, and x is shown in figure 1. The 

velocity v as shown in equation (10) is obtained for this and each 
subsequent case from the nomogram in figure 2 . 

Rectilinear flow between coaxial cylinders is shown in sketch (c). 



Diverging flow Converging flow 


(c) 

In 1913 a paper by Langmuir (ref. 3) was published that contained an 
approximate solution for this type of flow. The current density and 
current per unit length may be written 

j = 5. 467X10" 8 A -1 / 2 (cp o - cp) 3 / 2 r“ 2 p" 2 (27) 

y = 2jtr j = 3. 43X10" 7 A _1 / 2 (cp - cp) 3 / 2 (r|3 2 ) _1 (28) 

where cp is the potential at r Q and p is a nondimensional function 
of r expressed as a power series in ln(r/r Q ). 


In 1923 an improved solution hy Langmuir and Blodgett (ref. 4) was 
published. Two series representations for P were given. One of these 
most useful for electrode design is: 
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The Langmuir function 3 2 is plotted as a function of r/r Q in fig- 
ure 3 and is tabulated in table II. Note that 3 2 for r/r G >1 is 

for divergent flow (from the inner to the outer cylinder), and 3 2 for 
r/r Q < 1 is for convergent flow. The numbers in table II were computed 

from equation (29) using. all 14 coefficients given by Langmuir in refer- 
ence 4. These numbers agree well with those shown by Langmuir. Nomo- 
grams relating j/Z, (cp Q - cp), A, r, r/r Q , and v are shown in figures 
2, 4, and 5. 


Rectilinear flow between concentric spheres . - Sketch C also depicts 
rectilinear flow between concentric spheres. The solution for this case, 
due to Langmuir and Blodgett (ref. 5), was published in 1924. The cur- 
rent density and total current may be written 




J = 


= 5.467X10" 8 A -1 / 2 (cp o - cp) 3 / 2 r _2 a' 2 (30) 

4jtr 2 j = 6. 87x10" 7 A _1 / 2 (cp 0 - cp) 3 / 2 a" 2 (31) 


where a is a nondimensional function of r expressed as the power 
series 


a = pin ^-) - 0.3 (in ~^j Z + 0.075 (in ^-) 3 - 0.014 (in ~) 



The Langmuir function a 2 is plotted as a function of r/r Q in fig- 
ure 6 and is tabulated in table III. As in the cylindrical flow case, 
a 2 for r/r Q > 1 corresponds to diverging flow and a 2 for r/r Q < 1 

to converging flow. The numbers in table III were computed from equa- 
tion (32) and agree well with Langmuir r s data. 


Nomograms relating J, (cp 0 - cp), A, r/r Q , and v are shown in 
figures 2, 7, and 8. 
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Circular flow . - Circular flow is depicted in sketch (d). 



(a) 

The solution to this case is given by Meltzer in reference 11. The cur- 
rent density may he written 

j = 1.23X10" 7 A -1 / 2 (cp o - cp) 5 / 2 (r sin —j (33) 

where cp is the potential at (r, e) and j is the current density in 
the current sheet r. Figure 9 is a nomogram relating j, (cp Q - cp). A, 

r, and 9. The velocity can he obtained from figure 2. This flow has 
several peculiarities that are not evident from equation (33). Since 
these are discussed in reference 11, they will only he mentioned here. 
First of all, the ions accelerate through an angle of 60° and then de- 
celerate, reaching zero velocity at 120°. Since the equipotentials do 
not lie along a radial plane except at 0=0° and Q = 120°, the ve- 
locity varies along any radial plane in the flow. Finally, along any 

radial plane the current density j is proportional to l/r 5 . This 
poses a very difficult problem in ionizer design. 

If a beam is considered that is hounded hy an inner trajectory of 
radius a and an outer trajectory of radius h, the current per unit 
length may he written 

j = 3. 08x10" 8 A “ l/2 (^)(l - ^i)(cp 0 “ 9 ) 3/2 (^n 2 (34) 

This equation is used in the electrode design section. 
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Hyperbolic flow . - Hyperbolic flow is shown in the following sketch: 


y 



The solution to this case was obtained by Meltzer (ref. 7) and later by 
Rosenblatt (ref. 12). The space-charge-limited current density may be 
written: 

j = 4. 92x10” 7 A" 1//2 (cp 0 - cp) 3//2 r -2 (35) 

where cp is the potential at radius r. This equation is expressed in 
terms of a nomogram in figure 10. It is interesting to note that J 
along any given trajectory is a function only of the radial component of 
the trajectory. If the bounding trajectory is defined by the equation 
y = ( A/x), the current per unit length for one quadrant may be written 

j = 4. 92x10” 7 A _1 / 2 A(cp o - cp) 3 / 2 ^x 2 + (36) 

where Z is measured normal to the (x,y) plane. 
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Flow between two inclined planar electrodes , - The configuration of 
two inclined planar electrodes presents a very interesting problem of 
space-charge flow applicable to ion engine design. 

Walker (ref. 9) solved this problem by first assuming negligible 
space charge. In this case the potential distribution between two in- 
clined planar electrodes is known and can be written: 


By substituting (37) into Hamilton- Jacobi equation (9) in polar coordi- 
nates, equation (9) then becomes 


Because the right side of equation (38) is independent of r, so must be 
also the left side of (38). Therefore, Walker (ref. 9) sought a solution 
of the action function in the form 


By substituting (39) into (38) the Hamilton- Jacobi equation becomes 


(37) 


a 



(38) 


W = rz(0) 


(39) 



(40) 


where z = ^ and z ' = ^ Equation (40) was solved by Walker in 

the form of a power series. 


where z 


The components of velocity are then given by 

3 W 




Hence, 
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and after integration the trajectory* equation is obtained: 



where (r Q , 0) is the initial point of a particular trajectory, and (r, 0 ) 

are the coordinates of the remainder of that particular trajectory 
(sketch (f)). 



The angle between the trajectory and the normal to equipotential is 
given by 


tan cd = — (42) 
r d0 z T v 

In the case of space-charge flow between two inclined planar electrodes. 
Walker assumed an action function in the form 

w = rg(e) (43) 

By substituting (43) into Hamilton- Jacobi equation (9), equation (9) be- 
comes 


+ g ' 2 = 



$ 


.2 


(44) 


only; and. 
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Nov the left side of equation (44) is a function of 
hence, $ must he a function of 0 only. 


By expressing Poisson's equation (l) in polar coordinates, it 
follows that 



The components of vector current density are 

€ 0 $"g 

2 

r 

and 

G o*V 
3 e ~ r 2 


(45) 


By substituting these components into the equation of conservation of 
charge (eq. (3)) 


V 


A 


Sr 



i 

+ r = 


0 


Walker obtained 


± (®"g«) - $"g = 0 (46) 

By eliminating $ from equations (44) and (46) an ordinary fourth- 
order differential equation in g is obtained in the form 


+ g 


* ■) 



= 0 


(47) 


Any solution of this equation defines a possible motion. In a fashion 
similar to that in the space-charge-free case the trajectory equation 
is given by 


r 

— = exp 
r^ 



( 48 ) 
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and the angle between the trajectory and the normal to the equipotential 
is given by 


tan 


dr _ _g_ 
r d0 g 1 


(49) 


Ivey (ref. 10) made a detailed calculation of this problem and his re- 
sults of equations (4l), (42), (48), and (49) are presented in figures 
11 and 12. 


The space-charge-limited current density of the emitter as calculated 
in reference 10 is 


j 



/jkA 1 / 2 $3/2 

V m/ r 2 


F(a) 


(50) 


where F(a) could be defined as a perveance function depending on the 
magnitude of the angle a shown in sketch (f). At any place in the 
flow the current density may be written 


. jo^) 

J " r / r o 

Substitution into equation (50) yields 


j 



(51) 


where $ is the potential difference "between 6 = 0 and 6, and the 
range of 6 is 0 < 6 < a. A nomogram of equation (51) is shown in 

figure 13 where = 5. 46 7X10“ ^ A“^/^ for singly charged ions, 

y u f m 

The normalized potential distribution for space-charge-limited flow as 
a function of normalized angle 6 / a is shown in figure 14 for several 
values of a. It is interesting to note that for a = 0 the potential 
distribution is that for two parallel planes and, further, that the po- 
tential distribution for a < jt/2 is not greatly different from that 
for a = 0. 


Flow between two coaxial right circular conical electrodes whose 
vertexes coincide , - In this configuration it is also assumed, as in the 
case of the flow between two inclined planar electrodes, that the end 
effects can be taken care of in all directions. Walker has outlined the 
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procedure for the solution of this problem for the space-charge flow 
(see sketch (g)). The authors of this paper assumed an action function 
in the form 


W = r • s(e) (52) 



Kino and Harker (ref. 13) in their analysis used a similar form of action 
function W = r n • s(0) and investigated cases for n < 1. Some of the 
cases investigated in their paper show beam trajectories that pass through 
potential mi ni mums (accelerate-decelerate system). In the case of 
W = r • s(0), investigated herein, the beam trajectories are continuously 
accelerated (accelerate system). 

By substituting (52) into Hamilton- Jacobi equation (9) this equa- 
tion becomes 


s 2 + s' 2 



(53) 


As in the case of the flow between two inclined planes, it was deduced 
that cp has to be a function of 6 only. 
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By expressing Poisson's equation (l) in spherical polar coordinates, 
it follows that ’ 


= Sfe-0 + C P") 


(54) 


The components of the vector current density are 


^ + *") 


and 


e o s 


^6 r 2 Vtan 9 


(tin? + *') 


(55) 


The conservation of current equation (3) in spherical polar coordinates 
is 


V • j 


* /2J r 

]= l~ 




Je 


dr " l " r tan 6 + r d0 


i 




o 


(56) 


Substituting equation (55) into (56) gives then the condition for which 
div j is zero, namely, when 


cot 0(2cp"s' + cp’s") + cp'"s' + cp"s" - cp's' = 0 


(57) 


By eliminating cp from equations (53) and (57) a fourth-order nonlinear 
ordinary differential equation in s is obtained in the form 


s"" = - cot 6 


^(s"' 




+ s(l - ~T~j ~ 3s" (58) 


Equation (58) was solved numerically for s(0) and s'(0), which were 
then used for determination of ion trajectories. 


The components of velocity are 


dW 

v r = "Sr = 


s 
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and 



and 



and "by integration 


r 

— = exp 


o 



(59) 


The angle "between the trajectories and the normals to the equipotentials 
is given by 


The current density is given by 


To obtain more accurate values of j Q , the current density at the emitter, 
the current density j a at the collector was first calculated from 
equation (61). The current density at the emitter was then calculated 
from 


The solutions of ion trajectories and the angle between the trajectories 
and the normals to the equipotentials are presented in figures 15 and 16, 
respectively. A nomogram of equation (62) is given in figure 17. The 
normalized potential distribution as a function of normalized angle 
(e - a 1 )/(a 2 - c^) for several angles (a g - a x ) is shown in figure 18. 



(61) 



(62) 
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Partial Space -Charge Flow 

As early as 1920, Jaffe considered the plane diode under partial 
space- charge conditions where the current is less than space- charge 
limited, sometimes called temperature or flow-rate-limited condition, 
(ref. 14). Since that time there have been several papers on this sub- 
ject, two of which will he noted here. Brubaker (ref. 15) and Ivey 
(ref. 16) derived the equations governing the flow of electrons or ions 
under partial space- charge conditions independently at about the same 
time. Brubaker made use of normalized variables representing current 
density, voltage, distance, and emitter field, so that these variables 
range between zero and unity as the current density varies from zero to 
the space- charge-limited value. He defined the normalized variables as 

— Current density _ 

J ~~ Space-charge-limited current density j 

s 

— Voltage at distance X from emitter _ 

Voltage at collector, at distance X a from emitter ~ V a 

Y _ Distance from emitter _ X 
Emitter- collector spacing X a 

fa\/dV\ 

vU/x=o 

The previously mentioned assumption still holds, namely, zero velocity 
of the particles at the emitter. 

The potential distribution and electric field are fully defined by 
equations (63) and (64): 


— Gradient at emitter under influence of space charge 


T = 


Gradient at emitter in absence of space charge 


- 


16 Qv 3 / 2 - T^ 2 ) = 27(Vf 2 - Xf 3 ) 


(63) 


dV = 27T 5 - .32JX , (64 

to 27r 2 - 24JV 1 / 2 

The relation between the current density j and the_ele_ctric field at 
the emitter is given by equation (63). By setting V = X = 1 the con- 
dition at the collector is given by 


j -k * (i - f r )< 3r + 1)1/2 


(65) 
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Equation (65) is plotted in figure 19. It is interesting to see that 
the emitter field rises rapidly as the current is decreased from space- 
charge-limited values. Decreasing the current 10 percent from the space- 
charge -limited value causes the field to increase from zero to more than 
25 percent of its value in the absence of space charge. 


If_equation (63) is solved for_specific values of the electri£ 
field y and the current density j or, in other words, if y = y Q 

and j = j Q are held constant, equation (63) gives a relation between 
V and X. By substituting y 0 and j 0 into (63) and rearranging terms, 
equation (63) becomes 


X = 




1 + 


(«C 

\ f 




( 66 ) 


Equation (66) is plotted in figure 20. It is seen that the maximum de- 
pression of the space potential for space-charge-limited conditions 
(j 0 = l) occurs at x = 27/64 where the normalized slope is unity. By 

using the data obtained in figures 19 and 20 in equation (64) the elec- 
tric field is obtained as a function of position. This is presented in 
figure 21. The normalized gradient is equal to_unity in the plane 
x = 27/64 independent of the current density j . As Brubaker points 
out, this is only an approximation. However, the maximum deviation from 
unity is only 0.6 percent as the current density is varied from zero to 
space-charge-limited values. 

Ivey investigated the case of a cylindrical diode under partial 
space-charge conditions (ref. 16). In this case Poisson's equation can- 
not be integrated directly so a differential analyzer was used to obtain 
a numerical solution. Without a rigorous proof, Ivey then states that 
it appears that the emitter field characteristic as shown in figure 19 
is "universal" and applies also to other geometries, including those 
with external emitters. He points out also that this seems to be a di- 
rect consequence of Poisson's equation, and the only restrictions on the 
applicability of figure 19 are those of negligible initial velocities 
and an equipotential emitter (the same conditions imposed on the gener- 
alized space-charge law). 


THEORY OF PIERCE ELECTRODE DESIGN 

The problem of design of electrode systems that will constrain an 
accelerating ion beam to travel in a desired region has confronted sci- 
entists and engineers for many years. Before 1940 most electrode systems 
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were designed using assumptions of either zero space charge or paraxial 
flow. Since then much progress has heen made toward the development of 
general methods for the design of electrodes for high-charge-density ion 
teams. These methods can he divided into two categories: 

(1) Analog methods such as the resistance network and electrolytic 

tank 

(2) Analytic solution of the space-charge equations and of Laplace's 

equation in the region surrounding the ion team 

The most widely used analytic approach to accelerator design was 
proposed hy Pierce (ref. 17) in 1940. Although the Pierce method was 
originally intended for rectilinear flow, it has since heen generalized 
to include any type of curvilinear flow provided the boundary conditions 
are regular (ref. 18). When the solution to a particular type of space- 
charge flow (such as those described in the first section of this report) 
is found, it generally will fill all space. To form a beam it is neces- 
sary to take a segment of this flow and replace the missing portion with 
an electric field that will exactly match the conditions along the stream- 
lines bounding the ion beam. This is done by arranging suitably shaped 
electrodes around the ion beam. The electrode shapes are found by solv- 
ing Laplace's equation in the charge-free space surrounding the ion beam 
subject to the conditions of continuity of potential and normal gradient 
across the beam edge. 

When the space-charge flow can be expressed in terms of some type 
of cylindrical coordinate system, as is frequently the case, application 
of complex variable theory becomes quite useful. This is demonstrated 
in sketch (h). Curve C represents the boundary of a general curvilinear 
ion beam that extends to infinity normal to the paper. If an analytic 
function can be found which transforms C into the u-axis, then the po- 
tential distribution along C and the gradient normal to C will be trans- 
formed as shown. 
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Since harmonic functions remain harmonic under conformal transformation, 
a solution of Laplace's equation for the upper half of the W-plane can 
he converted to the desired solution in the Z-plane by application of the 
inverse transformation. A general solution of Laplace's equation in the 
upper-half plane subject to Cauchy boundary conditions is given in refer- 
ence 19: 


$(u, v) 


= Re$Q(W) + Im 





(67) 


•where 

Re real part 

Im imaginary part 

$ potential distribution along 

o 

K potential gradient normal to 
o 

In the event that N c cannot be integrated explicitly as an ana- 
lytic function of W, another technique described by Lomax (ref. 18) 
may be used to obtain $(u,v). 

It can be seen from equation (67) that $(u, v) can be written as 
the real part of a complex function: 



pW 

A = 0 + il = 3> C (W) - i / Nc(O d £ 


( 68 ) 


Hence, 


dA 

dW 


d$ c 

dF 


iN c (W) 


(69) 


By application of the Cauchy-Riemann conditions 


dA _ So _ 
dW "§u 'Sv 


"u 


+ iE y 


where E u and E y are the electric field components in the u and v 
directions transformed from the Z-plane by the mapping function W=f(Z). 
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Specifically, E u corresponds to the electric field tangent to the he am 
edge and E y corresponds to the field normal to the beam edge. 

Since E u and E v are derivable from an analytic function, they 

are differentiable; and the differential equation of a line of electric 
flux can be -written: 


dv 

du 



E u + ° 


(70) 


The equipotentials are orthogonal to the lines of force and are there- 
fore determined by the equation 


dv _ 
du E v 


E v ^ 0 


(71) 


The condition E v ^ 0 is always satisfied in cases of interest since 
E y = -Nq(u). By taking the real and imaginary parts of equation (69), 
equation (71) can be -written: 


dv 

du 




(72) 


This is the differential equation specifying the equipotential surfaces. 

While this technique is applicable to a large number of geometries, 
there are some configurations of interest that cannot be handled in this 
-way. An example of one of these is rectilinear flow with circular cross 
section. The solution of this particular problem will be discussed 
later. 


Electrodes for Rectilinear Flow (Nq = 0) 

Rectilinear slab beam . - This is an application of the solution of 
flow between infinite parallel planes to the formation of a finite slab 
ion beam as shown in sketch (i). 
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Mathematical model 

(i) 


From equation (25) the potential distribution along the beam edge is 


$> 


C 



(73) 


C-l = 6.942X10 4 A 1 / 3 ,) 2 / 3 


There is no need for conformal mapping . since the beam edge already co- 
incides with the real axis. The boundary conditions to be satisfied are 
shown in sketch (i). The potential gradient normal to the beam edge Nq 

is zero. Therefore, equation (67) has the particularly simple form: 

4/3 

$ = ReO c (Z) = C^ReZ ' 

= C-j_(x 2 + y 2 ) 2 ^ 3 cos^| 'tan -1 ^ (74) 

Equation (74) is developed in dimensional form for clarity. It is much 
more convenient, however, to nondimensionalize from the start. This may 

be done by letting X = x/Z, Y = y/Z, and X = ^^/(C^Z 4 / 3 ) where Z is 

the unit of length. Then equation (73) becomes 


( 75 ) 
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Sketch (i) remains the same except that now the boundary conditions are: 



(75) 



y=o 


= 0 


and equation (74) becomes 


X = Re(X + iY) 4 / 3 = (X 2 + Y 2 ) 2 ^ 3 cos(| tan" 1 |) (76) 

Equipotentials defined "by equation (76) are plotted in figure 22 and 
tabulated in table IV. Note that the X = 0 electrode makes an angle 
of 67.5° with the beam edge. This is called the Pierce angle , and it 
occurs in every geometry because, as the ion emitter is approached 
closer and closer, the potential distribution approaches that of a plane 
diode. The function X 4 / 3 is multiple -valued and hence can give rise 
to solutions that do not satisfy the boundary conditions. While it is a 
simple matter to choose the proper solution in this case, it is not so 
simple in some of the cases that follow where the complex potential func- 
tion is a power series. The coordinates listed in table IV agree with 
the curves shown by Pierce in reference 17. 

Beams of circular cross section . - The potential distribution has 
the same form as equation (73), but in this case polar coordinates are 
used and x is replaced by z: 

Oq = 

C ± = 6.942X10 4 A 1 / 3 ^ 2 / 3 (77) 

In 1955 a paper by Daykin (ref. 20) appeared in which a solution to this 
problem was described. A sketch of the problem is shown as follows: 



tion 

U) 
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Equation (77) can "be nondimensionalized "by writing 


*C = 


°€ 



(78) 


where a is the radius of the ion beam. Daykin showed that the poten- 
tial distribution near the beam edge can be expressed by an infinite 
series of the form 



The functions F 2 and F 4 are plotted in figure 23. Because of its 

rapid convergence near the beam edge, this series may be truncated at 
three terms without introducing serious error. In order to exbend the 
solution away from the beam edge the asymptotic forms of the equipoten- 
tials are computed. For this purpose it is convenient to define a new 
coordinate system in terms of r and z. This is shown in sketch (k). 


r 



1/2 


The equipotentials are then written in terms of Legendre functions P 4 / 3 
of order 4/3 as 

/■d\ 4 /3 

X = (-) ^4/3^ cos (80) 

Equation (80) is used to obtain electrode shapes far from the beam edge. 
Equipotentials from equations (79) and (80) are plotted in figure 24(a) 
and tabulated in tables V and VI. The transition region is shown by 
dotted lines. 

For a hollow ion beam another set of electrodes is required inside 
the beam. Daykin derived an expression for these equipotentials for 
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small r/b where h is the inner radius of the beam. This equation may- 
be written 



As in the previous case it is necessary to fair in part of the equipo- 
tentials. Some of these are shown in figure 24(b) and tabulated in 
table VII. The dotted portion joins the equipotentials from equa- 
tion (81) to the proper intercept on the inner beam edge defined by 
equation (78). Until recently this was the only analytic solution of 
the Pierce electrode problem for axisymmetric ion beams. Harker (ref. 21) 
has recently developed a very general method whereby any axisymmetric 
electrode design problem can be solved with the aid of a digital computer. 

Cylindrical flow . - The application of converging and diverging flow 
between coaxial cylinders is considered in this section. The two con- 
figurations are shown in sketches (l) and (m). 

Converging flow 




Electrode shapes 


( 1 ) 
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Diverging flow 




(m) 


Since the included angle of the wedge is a, the current per unit length 
from equation (28) will he multiplied by the factor a/2rt. The potential 
distribution becomes 


C 2 


0 C = C 2 (r|3 2 ) 2 / 3 

2 /3 

6. 9 42x10 4 A 1 / 3 ^-) 


(82) 


This is nondimensionalized by letting R = r/r Q and X^ = ^/(Cgr^/ 3 ) 
The nondimensional potential distribution becomes 


2/3 


X = (R(3 2 ) 


(83) 
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From equation (29) it follows that 


(rp^) = 


—I 


V B n (ln B) n 
n=l 


and 


X = 


-00 

X) B » <in E > n 


Ln=l 


4/3 


(84) 


Since the potential gradient Nq normal to the ion beam is zero, the 
expression for the potential outside of the beam is, from equation (67), 


where 



(85) 


Z 


r 


o 



10 


Re 1 ® = R(cos 0 + i sin 6 ) 


Equation (85) can be rewritten 


X = (P 2 + Q 2 ) 2// ° cos[| tan" 1 (86) 

where 

P = [in R + 0.l(ln 2 R - 0 2 ) + 0.0167(ln R 3 - 30 2 In R) 

+ 0.00242(ln 4 R - 60 2 In R + 0 4 ) + . . .] 

and 

Q = [e + 0.2 e In R + 0.0167(30 ln 2 R - 0 3 ) 

+ 0.00242(40 lipR - 40 3 In R) + . . .] 
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Equation (86) was solved for several values of X on a digital computer. 
Equipotentials for both converging and diverging flow are shown in fig- 
ures 26 and 27 and are tabulated in tables VIII and IX, respectively. 
This case was also solved by Radley (ref. 22) for convergent flow, and 
the results shown in table VII agree with those shown in reference 22. 


Circular flow. - The space-charge-limited- flow equations for which 
the ion trajectories are concentric circular arcs are used here in ^ the 
design of a curved ion beam. This solution was obtained by Lomax in 
reference 18. The potential distribution obtained from equation (34) is: 


Electrodes for Curvilinear Flow (W c / 0) 



(87) 


where 



The potential gradient normal to the beam is 



( 88 ) 


A sketch of the electrode structure is shown as follows: 



(n) 
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From the sketch it is seen that the problem can be divided in tvo parts: 
r > b and r < a. Sketches for these two problems together with the 
corresponding boundary conditions are shown as follows: (The shaded 

region is the region of interest in each case. ) 



(o) 


In this case it is necessary to map the beam edge onto the real axis. 
Lomax shows that a mapping function that will do this is (for r = b) 

Z = be lW (89) 

The boundary conditions become, on v = 0, 



The formal solution of Laplace's equation in the upper half of the W- 
plane is given by equation (67), which becomes 
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<J>(u, v) 



(90) 


Explicit integration of the second, term is not readily carried out; 
therefore, the differential equation representing the equipotentials 
must he solved. By combining equations (69), (72), and (90), the dif- 
ferential equation becomes 


dv 

du 



(91) 


From the transformation equation (89) it follows that 



e 


-v 


and 

6 = u 

Equation (91) becomes in the Z-plane 



Equation (92) was integrated using the Runge-Kutta method, and some 
typical equipotentials are shown for a/b =0.7 in figure 27 and are 
tabulated in tables X and XI. The inside equipotentials were found by 
simply replacing b by a in equations (89) and (92). In figure 27 

X b = $(b 2 /C 3 ) and X a = $(a 2 /C 3 ). These results agree with those ob- 
tained "by Lomax in reference 18. 

Electrodes for hyperbolic flow . Sketch (p) shows a segment of 
hyperbolic flow bounded by the curves y = l/x and 10/x. This solu- 
tion was obtained by Rosenblatt in reference 12. 
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The potential distribution along the beam edge as given in the space- 
charge- flow section is 


0 


C 
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where 


C 4 = 1.60X10 4 A 1 / 3 ^ ^ 


and 


A = 1,10 

The solution of Laplace's equation outside the ion beam in the upper 
part of the first quadrant (i.e., outside the boundary A = 10) will he 
obtained first. It is shown in reference 12 that a transformation that 
will map the beam boundary onto the real axis is: 


W = i (Z^ - 2iA) = u + iv ^ 


U = | (X 2 - y 2 ) 




(93) 


v 


xy - A 
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A sketch of the transformed boundary is shown as follows: 


v 



(q) 

The transformed "boundary conditions are also shown. Equation (67) may 
now "be used: 

/»W 

cp(u,v) = Re£-2A(W 2 + A 2 )^] + Im / |^-2C 4 A(£ 2 + A 2 )’ 1 ^^ 

Jo 


0 

After integration and separation of the real and imaginary parts the 
equipotentials in the W-plane are expressed by: 


(94) 


®(u,v) = -2C 4 |-^r {[(u 2 + A 2 - v 2 ) 2 + 4u 2 v 2 J 1//2 + (u 2 + A 2 - v 2 )| 


1/2 


v + 


1 / 2 ' 


+ A tan 


-1 


j[( U 2 + A 2 - y 2 ) 2 + 4u 2 t 2 ] 1/2 - (u 2 + A 2 - v 2 )} 
u + -i={[(u 2 + A 2 - v 2 ) 2 + 4 u 2 v 2 J 1/2 + (u 2 + A 2 - v 2 )} 


( 95 ) 
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This equation may he nondimensionalized hy dividing equation (95) hy 
-C^A. Substitution of the values of u and v from equation (93) yields 

the expression for the electrode shapes in the (x,y) plane. Some of 
these are plotted in figure 28. Equipotentials outside the beam edge 

A = 1 may be obtained by the mapping W = i (2iA - Z 2 ). The potential 

function is identical to equation (95) except for minus signs in front 
of the nonsquared u and v in the arc tangent term. Some electrode 
shapes for this boundary are also shown in figure 28. 

Since there is no zero velocity surface in this configuration, it 
may not have any practical application. It may possibly be used to de- 
flect an ion beam, however. 


Limitations of the Pierce Technique 

The departure of analytically derived ion accelerators from physical 
fact lies mainly in two areas. The first of these is that practical re- 
quirements dictate modifications in the calculated electrode shapes. 

These modifications manifest themselves in the curtailment of electrode 
size and the insertion of apertures in the electrodes to pass the ion 
beam. The inherent instability in the solution of elliptic differential 
equations with Cauchy boundary conditions on open boundaries works in 
reverse to minimize the error caused by curtailment or modification of 
electrode shapes away from the boundary. This means that, while small 
variations of potential along the boundary will cause increasingly larger 
variations the farther away from the boundary, conversely, large varia- 
tions far away from the boundary cause only small variations along the 
boundary. This is readily observed in an electrolytic tank where a 
tremendously wide variety of electrode shapes will give the proper po- 
tential distribution along the beam edge to within the accuracy of the 
measuring equipment. It is difficult to give exact numbers for the 
magnitude of this effect other than to say that the error propagates 

away from the boundary approximately as e Ky , where y is the distance 
measured from the boundary and K is a constant. 

Aperture effect . - Insertion of apertures in the electrodes has two 
deleterious effects on ion accelerator performance. The first is that 
the aperture acts as an electrostatic lens which will tend to diverge an 
accelerating ion beam. If multiple electrodes are being used, this re- 
sults in impingement of ions on the electrodes. Even if only one accel- 
erating electrode is used, the formation of a diverging ion beam will 
result and may cause some impingement on this electrode. It is shown in 
reference 23 that the focal length of an aperture is 

, __ ^2g0 


(96) 
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where g = 1 for slits, g = 2 for circular apertures, and is the 

negative potential gradient on the upstream (or source ) side of the 
aperture. The meaning of this is demonstrated in sketch (r). 



The second effect is to reduce the ion current by increasing the 
effective accelerating distance. Depending on the ratio of aperture 
width to accelerating distance, experimentally observed current densities 
can he a small fraction of the theoretically expected values. 

Emitting surface irregularities . - The second area of departure of 
the Pierce theory from physical fact lies in the assumptions regarding 
initial conditions at the ion emitter. The usually assumed condition of 
zero initial velocity along a perfectly smooth equipotential surface is 
not met in practice, and it will he shown that considerable error can 
arise as a result of this. Rather than treat these effects separately, 
the nature of the ion trajectories in the neighborhood of a rough emitter 
will be investigated, and dispersion due to transverse electric fields 
will be compared directly with thermal spreading. In order to do this 
it is first necessary to examine the potential distribution near a sur- 
face irregularity. A sketch of the mathematical model is shown as 
follows : 
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The surface irregularity is a cylindrical hump of radius Z, and it 
is assumed that at 100 units above the hump the equipotentials are es- 
sentially plane. An approximate expression for a flux line as computed 
in appendix B is (u/z) = (x/z){l + l/[(x/z) 2 + (y/Z) 2 ]f. The bounding 
flux line, which originates at the point x = Z, y = 0, is (u/z) = 2. 
This line intersects the line y = 100Z at x ^ 21. Hence, the flux 
line spreading due to the electrostatic field of the hump is Ax « Z. 

It is shown in appendix B that thermal spreading Axip at y = 100Z, 

which is the product of the most probable therma l velocity and the ion 
transit time, can be written Axip = 200Z-y/0. l/cp 0 where the effective 

thermal potential is assumed to be 0.1 volt. The ratio of electrostatic 
to thermal spreading is 


Ax _ / ^o 

^T \ 4xl0 3 


(97) 


The macroscopic electric field strength E in an ion rocket is defined 
as the applied accelerating voltage divided by the accelerator spacing. 
For emission-limited flow E usually has a value of between 10^ and 10 7 
volts per meter. From sketch (s) it can be seen that the voltage across 
the model cp Q can be written cp Q = E(100Z), and the spreading ratio be- 
comes 


Ax 

Axip 



(98) 


Equation (98) is shown plotted in figure 29. It is seen from figure 29 
that spreading due to the electrostatic field of a surface irregularity 
can be of the same order of magnitude as thermal spreading. It should 
be noted that this analysis is based on zero space charge and therefore 
is only approximate. 


Analysis of accidental misalinement of parallel planar electrodes . - 
As pointed out by Ivey (ref. 10), an accidental misalinement of parallel 
planar electrodes can be analyzed as an application of the motion be- 
tween two inclined planar electrodes. Consider the electrode arrangement 
as shown in sketch (t): 
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One electrode is tilted about its center by an angle a from the paral- 
lel condition and 


sin a 

The current per unit length when a = 0 may be written 


J cc=0 _ 4 /^qN 1 / 2 <D 3 / 2 w 

l 9 e o[ m ) s 2 


(99) 


and thus the ratio J/J x=0 , which gives the effect of the tilt, is 


J 

J a=0 



( 100 ) 


Values from equation (100) were calculated by Ivey for small angles of 
a and are presented in figure 30. For a rectangular emitter with an 
aspect ratio of w/s = 1.24, the current per unit length is insensitive 
to the angle of tilt about a line passing through its center and parallel 
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to its side. Ivey also showed the independence of angular misalinement 
in an arbitrary direction by the use of a disk- shaped emitter. The total 
current for a = 0 in this case is 




1/2 j.3/2,.,2 


4s 


( 101 ) 


where D is now the diameter of the disk, and the effect of electrode 
tilt becomes 



Values from equation (102) were also calculated by Ivey and are presented 
in figure 31. It is seen that the variation of current for the aspect 
ratio D/s = 1.42 is zero. For this optimum geometry (aspect ratio 
D/s = 1.42) equation (101) becomes 


1/2 


J^ 0 = 6.233X10-1 2 ^ 0 3/2 


(103) 


and the behavior of J/J^q for "the above optimum geometry is shown in 

figure 31. It is seen that even a large angle of tilt does not change 
the total current appreciably. This, of course, does not mean that the 
ion optics will not be affected, and in an ion rocket the ion impinge- 
ment on the electrodes must be reduced to nearly zero for long-duration 
operation. Therefore, electrode alinement is still quite critical. 


CONCLUDING REMARKS 

Three methods for the analytic solution of the space- charge -f low 
problem appear in the literature, in which certain properties of the 
flow are prescribed. These are (l) initial and final equipotentials, 

(2) velocity and acceleration, and (3) trajectories and action function. 
None of these methods is generally applicable at present. Of the three 
methods, the trajectory-action function method appears the most versatile 
for the design of electrostatic rocket engine accelerators. 

The existing solutions of space- charge flow have been described, 
and the parameters pertinent to the design of electrostatic rocket engine 
accelerators have been given for the following geometries: (l) paraxial 
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flow, (2) flow between coaxial cylinders, (3) flow between concentric 
spheres, (4) circular flow, (5) hyperbolic flow, (6) flow between in- 
clined planes, and (7) flow between coaxial cones. 

Electrode shapes determined by the method of Pierce have been given 
for the following geometries: (l) rectilinear flow with rectangular 

cross section, (2) paraxial flow with circular cross section, (3) con- 
vergent flow between coaxial cylinders, (4) divergent flow between co- 
axial cylinders, (5) circular flow, and (6) hyperbolic flow. 

Limitations of the space-charge-flow and electrode design theories 
have been discussed. One is the aperture effect and another is the ef- 
fect of irregularities on the ion emitting surface. Surface irregular- 
ities have been demonstrated to cause lateral ion velocity components of 
the same order of magnitude as thermal ion velocities. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 20, 1962 
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APPENDIX A 
SYMBOLS 

(inks units used throughout. ) 

A 

molecular weight/ionic charge (ionic charge may he greater 
than one in multiply charged ions) 

a 

inner radius of circular ion he am; outer radius of ion 
heam with circular cross section 

B 

(h 2 h 3 /h 1 )v 2 ( 1 /h|) 

h 

outer radius of circular ion heam; inner radius of ion 
heam with circular cross section 

C 

h 2 h 3 /h^; also ion heam edge 

C 1' C 2> C 3' C 4 

constants identified in text 

D 

diameter of disk emitter 

d 

S/sin a, for electrode tilt analysis 

E 

electric field strength, volts/meter 

e 

unit electronic charge, 1.602x10“-^ coulomb 

F 2; F 4 

dimensionless functions of r/a defined in ref. 16 

F(a) 

function defined in eq. (50) 

g(0) 

function defined in eq. (43) 


gravitational conversion factor, 9.806 meters/sec^ 

H(a) 

function defined in eq. (100) 

h l> h 2> h 3 

scale factors in curvilinear coordinates, 
ds = h^ d£ + h 2 d*rj + h^ d£ 

Im 

imaginary part 

J 

total current, amp 
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J/Z 

j 

k(ti,£) 

l 

m 

N c 

q 

R 

Re 


S 

s(e) 

ds 


u, v 

V 


V v y' v z 


v r> v 6 


v 

W 


w 

X,Y 


current per unit length, amp/meter 
current density, amp/meter 2 

(h 2 h 3 /e 0 )jW 2 q 

accelerator height, meters; also unit of length 
particle mass, kg 

potential gradient normal to boundary C 
charge, coulombs 

r / r o> r / r a* r / r t 
real part 

radius as in spherical and cylindrical polar coordinates 

emitter radius; also initial radius of trajectory in flow 
between inclined planes; also inner radius in circular 
flow engine 

accelerating distance for electrode tilt analysis 
function defined in eq. (52) 
differential displacement along a trajectory 
coordinates in W-plane (W = u + iv) 
ion velocity 

velocity components in Cartesian coordinates 
velocity components in polar coordinates 

acceleration vector 

most probable thermal velocity 

action function, meter 3 / sec; also complex number u + iv 
tilt analysis 

nondimensional Cartesian coordinates 
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x,y, z 
y(e) 
z 
z 

a 


a-pCCg 

P 

r 


r 


€ 


O 


e 

A 

A 

I 

p 

0 

$c 

cp 

*0 

X 


Cartesian coordinates 
function defined in eq. (39) 
x + iy 
x - iy 

Langmuir function for spheres; angle between inclined 
planes; arbitrary fixed angle 

conical half- angles 

Langmuir function for cylinders 

h^$ 

angle between ion trajectory and flux line in flow between 
coaxial cones 

-12 

permittivity of free space, 8.855x10 
coulomb/( newton )(meter^ ) 

curvilinear coordinate defining trajectories 

angular displacement, deg 

complex potential function, 0 + i¥ 

hyperbolic trajectory parameter 

curvilinear coordinate defining action function 

charge density, coulomb/meter 3 

cp o - cp 

potential distribution along beam boundary 

angular displacement in spherical coordinates; also poten- 
tial, volts 

emitter potential, volts 
nondimensional potential distribution 


Y 

CD 

Superscripts: 
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flux line 

angle "between ion trajectory and flux line for flow between 
coaxial cones 


T 11 Til 


If IT 


differentiation with respect to variable 


0 
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APPENDIX B 


APPROXIMATE ION BEAM DISPERSION DUE TO SURFACE ROUGHNESS 

In this appendix the electric field in the neighborhood of a cylin- 
drical hump will be obtained. Following that the transit time of an ion 
through the diode will be calculated, and an expression for thermal 
spreading, which is the product of the most probable thermal velocity 
and the ion transit time, will be shown. 

iy iv 



(u) 


The surface is an infinite plane containing a cylindrical hump of radius 
l. It is assumed that at a distance y = 100Z the equipotential sur- 
faces are plane. By the conformal transformation w = z + l/z the sys- 
tem is converted into two plane electrodes separated by a distance of 
99.99 l. This is shown on the right side of sketch (u). The space- 
charge-free potential distribution in the W-plane is 

<p(u,v) = cp o (l - 99> 99 ( B1 ) 


The electric field is obtained by differentiation of equation (Bl) with 
respect to v: 


(B2) 
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Returning to the right side of sketch (u) the flux lines are given 
by u /l = constant. From the transformation equation the expression for 
u /l is 


u 

7 



(B3) 


This equation is used in the computation of spreading due to the electric 
field of a surface irregularity. 

In the computation of transit time through the diode in sketch (u) 
(left side), it is assumed that the electric field is constant; that is. 


E = 


<Po 

1007 


(B4) 


The equation of motion of an ion is 


d 2 y _ „ W 0 


m — p = qE = 
dt 


(B5) 


Integration of equation (B5) subject to the initial conditions of 

•^ = 0 at t = 0 
dt 


and 


yields 


y = 0 at t = 0 


Wo ,2 


7 ~ 200m7 

At y = 1007 the transit time is 


t = 1007 




(B6 ) 


(B7 ) 
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From reference 24 it is seen that the most probable thermal velocity is 


v = 



(B8) 


where cp is the mean thermal potential (usually taken to he about 0.1 
volt). Thermal spreading is then the product of equations (B7) and (B8): 


Axip = 2002 




(B9) 
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APPENDIX C 

RESCALING OF NOMOGRAMS 

The nomograms shown in this report may he easily rescaled to allow 
computations "beyond the range of one or more of the variables. Consider, 
for example, the nomogram in figure 1, which represents the space-charge- 
flow equation for a plane diode. The equation is 

j = 5. 467x10“® A _1 / 2 (cp 0 - cp) 3 / 2 x -2 (26) 

The ranges of the various scales are 

10 < A < 10 3 

2 X 10 < J < 10 3 

10 “ 4 < x 5 10-2 

!0 2 < (cp Q - cp) < 10 4 

Suppose that it is desired to calculate the current density for heavier 
particles, say up to 10 5 amu. This means that the new range of the A 
scale will he from 10 3 to 10 5 amu. It is seen from equation (26) that, 
if (cp Q - cp) and x are held constant and if A is multiplied hy 10 2 , 

then J will he divided hy 10. That is, there exists a scale factor 
associated with A which, when multiplied hy the range of J, will give 
the new range. In this case, K A = 10" 1 and the new range of J will 

he 

2 < J < 10 2 

Usually it is necessary to rescale several variables simultaneously. 

This may he done in a similar manner. Suppose, for instance, that the 
desired ranges are to he changed in accordance with the following list: 

10 < A < 10 3 -»• 10 3 < A < 10 5 
10 3 -► ? < J < ? 


2 X 10 < J < 
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10" 4 < x < 10" 2 -» 10" 2 < x < 10° 

10 2 < (cp o - cp) < 10 4 -► 10 4 < (cp Q - cp) < 10 6 

Referring again to equation (26) it is seen that scaling factors may he 
assigned to each variable. For instance, it was shown that the scaling 

factor due to a hundredfold increase in A is lO*" 1 . The net scal- 

ing factor applied to J will he the products of the individual factors 
These are 

% = i 0 ' 1 

K x = 10-4 
K<p - 10 3 

and finally 

K J * Wq> - 10-2 

Therefore, the new range of the J scale will he 

2X10" 1 < J < 10 
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TABLE I. - SPACE-CHARGE -FLOW EQUATIONS FOR SEVERAL GEOMETRIES 


Type of flow 

i 

T1 

C 

^1 

hg 

h 3 

B 

C 

K(T1,0 

General equation (24) 

Rectilinear 

between 

plane 

electrodes 

X 

y = const. 

z = 
const. 

1 

1 

1 

0 

1 

6 0 T2q 

rVa A-k(h.c) 

dx^ 

Rectilinear 

between 

concentric 

cylindrical 

electrodes 

r = (x 2 + y 2) 1 / 2 

0 = cos-1 X 

z = 
const. 

1 

r 

1 

0 

r 

si J % 

«o » 2 i 

( r S?)] = 

Rectilinear 

between 

concentric 

spherical 

electrodes 

(x 2 + y 2 + Z 2 ) 1/2 

- -1 Z 

0 = COS p 

<P = 
const. 

1 

r 

r sin 0 

0 

r 2 sin 0 

r g (sln fl)J JUT 
£ o V2< ! 

1,1/2 (r 2 sin 0 - K(1U) 

Circular 

0 - cos -1 -p 

r = 

(X 2 + y 2 ) l/2 

z = 
const. 

r 

1 

1 

4r" 5 

r- 5 

J- iff 

e 0 T2q 

* V/2 [S + 4r ] - k <™> 

Hyperbolic 

|(* 2 - y 2 ) 

xy 

z = 
const. 

V 2 .... 

1/2 . 

1 

4 

4(? 2 + n 2 ) 

J Vm/2q 

4 ^ /2 [l (52+T ^ ) (tc) +r = 

U 2 + n 2 ) 1/2 

(? 2 + Tl 2 ) 1/2 

2e 0 U 2 + T1 2 ) 172 


VI = 0.7199X10" 4 VA. 


Identical to Meltzer's equation in ref. 10. 

*** 

This equation Is treated in ref. 11 and in a different fashion in ref. 6. 
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TABLE II. - LANGMUIR FUNCTION FOR CONCENTRIC CYLINDRICAL FLOW 


r / r o 

P 2 

r / r o 

3 2 

r / r o 

3 2 

r A 0 

3 2 

O.OIOOO 

1174.9000 

0.2500 

6.06010 

1.00 

0.00000 

4.0 

0.6671 

.01111 

1018.5000 

.2632 

5.37950 

1.01 

.00010 

4.2 

.6902 

.01250 

867 .1100 

.2778 

4.72980 

1.02 

.00039 

4.4 

.7115 

.01429 

721.4300 

.2941 

4.11260 

1.04 

.00149 

4.6 

.7313 

.01667 

582.1400 

.3125 

3.52930 

1.06 

.00324 

4.8 

.7496 

.02000 

450.2300 

.3333 

2.98140 

1.08 

.00557 

5.0 

.7666 

.02500 

327.0100 

.3448 

2.72140 

1.10 

.00842 

5.2 

.7825 

.03333 

214.4200 

.3571 

2.47080 

1.15 

.01747 

5.4 

.7973 

.05000 

115.6400 

.3704 

2.23010 

1.20 

.02875 

5.6 

.8111 

.05556 

97.9970 

.3846 

1.99950 

1.30 

.05589 

5.8 

.8241 

.06250 

81.2030 

.4000 

1.77920 

1.40 

.08672 

6.0 

.8362 

.07143 

65.3520 

.4167 

1.56970 

1.50 

.11934 

6.5 

.8635 

.08333 

50.5590 

.4348 

1.37120 

1.60 

.15250 

7.0 

.8870 

.10000 

36.9760 

.4545 

1.18400 

1.70 

.18540 

7.5 

.9074 

.10530 

33.7910 

.4762 

1.00860 

1.80 

.21770 

8.0 

.9253 

.11110 

30 . 6980 

.5000 

.84540 

1.90 

.24910 

8.5 

.9410 

.11760 

27.7010 

.5263 

.69470 

2.00 

.27930 

9.0 

.9548 

.12500 

24.8050 

.5556 

.55720 

2.10 

.30830 

9.5 

.9672 

.13330 

22.0150 

.5882 

.43320 

2.20 

.33610 

10.0 

.9782 

.14290 

19.3370 

.6250 

.32330 

2.30 

.36260 

12.0 

1.0122 

.15380 

16.7770 

.6667 

.22820 

2.40 

.38790 

14.0 

1.0352 

.16670 

14.3430 

.7143 

.14856 

2.50 

.41210 

16.0 

1.0513 

.17240 

13.4070 

.7692 

.08504 

2.60 

.43510 

18.0 

1.0630 

.17860 

12.4930 

.8333 

.03849 

2.70 

.45710 

20.0 

1.0715 

.18520 

11.6010 

.8696 

.02186 

2.80 

.47800 

30.0 

1.0908 

.19230 

10.7330 

.9091 

.00980 

2.90 

.49800 

40.0 

1.0946 

.20000 

9.8887 

.9259 

.00630 

3.00 

.51700 

50.0 

1.0936 

.20830 

9.0696 

.9434 

.00356 

3.20 

.55260 

60.0 

1.0910 

.21740 

8.2763 

.9615 

.00159 

3.40 

.58510 

70.0 

1.0878 

.22730 

7.5096 

.9804 

.00040 

3.60 

.61480 

80.0 

1 .0845 

.23810 

6.7705 

.9901 

1.0000 

.00010 

.00000 

3.80 

.64200 

90.0 

100.0 

1.0813 

1.0782 
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TABLE III. - LANGMUIR FUNCTION FOR CONCENTRIC SPHERICAL FLOW 


r / r o 

a ? 

r / r o 

a ? 

r / r o 

a 2 

r A 0 

a ,2 

O.OIOOO 

1144.000 

0.2500 

4.9680 

1.00 

0.0000 

4.2 

0.979 

.01111 

974.100 

.2632 

4.4290 

1.05 

.0023 

4.4 

1.022 

.01250 

813.700 

.2778 

3.9130 

1.10 

.0086 

4.6 

1.063 

.01429 

663.300 

.2941 

3.4210 

1.15 

.0180 

4.8 

1.103 

.01667 

523.600 

.3125 

2.9540 

1.20 

.0299 

5.0 

1.141 

.02000 

395.300 

.3333 

2.5120 

1.25 

.0437 

5.2 

1.178 

.02500 

279.600 

.3448 

2.3020 

1.30 

.0591 

5.4 

1.213 

.03333 

178.200 

.3571 

2.0980 

1.35 

.0756 

5.6 

1.247 

.05000 

93.240 

.3704 

1.9010 

1.40 

.0931 

5.8 

1.280 

.05556 

78.560 

.3846 

1.7120 

1.45 

.1114 

6.0 

1.311 

.06250 

64.740 

.4000 

1.5310 

1.50 

.1302 

6.5 

1.385 

.07143 

51.860 

.4167 

1.3580 

1.60 

.1688 

7.0 

1.453 

.08333 

39.980 

.4348 

1.1930 

1.70 

.2080 

7.5 

1.516 

.10000 

29.190 

.4545 

1.1036 

1.80 

.2480 

8.0 

1.575 

.10530 

26.680 

.4762 

.8880 

1.90 

.2870 

8.5 

1.630 

.11110 

24.250 

.5000 

.7500 

2.00 

.3260 

9.0 

1.682 

.11760 

21.890 

.5263 

.6210 

2.10 

.3640 

9.5 

1.731 

.12500 

19.620 

.5556 

.5020 

2.20 

.4020 

10.0 

1.777 

.13330 

17.440 

.5882 

.3940 

2.30 

.4380 

12.0 

1.938 

.14290 

15.350 

.6250 

.2968 

2.40 

.4740 

14.0 

2.073 

.15830 

13.350 

.6667 

.2118 

2.50 

.5090 

16.0 

2.189 

.16670 

11.460 

.6897 

.1740 

2.60 

.5430 

18.0 

2.289 

.17240 

10.730 

.7143 

.1396 

2.70 

.5760 

20.0 

2.378 

.17860 

10.010 

.7407 

.1084 

2.80 

.6080 

30.0 

2.713 

.18520 

9.315 

.7692 

.0809 

2.90 

.6390 

40.0 

2.944 

.19230 

8.636 

.8000 

.0571 

3.00 

.6690 

50.0 

3.120 

.20000 

7.976 

.8333 

.0372 

3.20 

.7270 

60.0 

3.261 

.20830 

7.334 

.8696 

.0213 

3.40 

.7830 

70.0 

3.380 

.21740 

6.712 

.9091 

.0096 

3.60 

.8360 

80.0 

3.482 

.22730 

6.109 

.9524 

.0024 

3.80 

.8860 

90.0 

3.572 

.23810 

5.528 

1.0000 

.0000 

4.00 

.9340 

100.0 

3.652 


TABLE IV. - ELECTRODE COORDINATES 
FOR RECTILINEAR SLAB ION BEAM 


Y 

X 

X 

X 

X = 0 

X = 0.4 

X = 1.0 

0 

0 

0 

1.0000 

.1 

.04141 

.50626 

1.0017 

.2 

.08283 

.51578 

1.0066 

.3 

.12425 

.53072 

1.0147 

.4 

.16567 

.55009 

1.0258 

.5 

.20708 

.57298 

1.0395 

. 6 

.24850 

.59864 

1.0558 

.7 

.28992 

.62648 

1.0743 

.8 

.33133 

.65607 

1.0947 

.9 

.37275 

.68708 

1.1169 

1.0 

.41416 

.71924 

1.1406 

2.0 

.82833 

1.07580 

1.4338 

3.0 

1.24250 

1.46010 

1.7793 

4.0 

1.65670 

1.85490 

2.1477 

5.0 

2.07080 

2.25520 

2.5285 

6.0 

2.48500 

2.65870 

2.9167 

7.0 

2.89920 

3.06440 

3.3100 

8.0 

3.31330 

3.47150 

3.7067 

9.0 

3.72750 

3.87970 

4.1061 

10.0 

4.14160 

4.28870 

4.5075 


TABLE V 


ELECTRODE COORDINATES FOR 


RECTILINEAR ION BEAM WITH CIRCULAR 
CROSS SECTION NEAR BEAM EDGE 


R 

z/a 


X = o 

X = 0.5 

X = 1.0 

X = 2.0 

1.00 

0 

0.593 

1.000 

1.683 

1.25 


.611 

1.009 

1.687 

1.50 


.649 

1.033 

1.702 

1.75 


.699 

1.072 

1.727 

2.00 


.751 

1.115 

1.756 

2.25 


.804 

1.161 

1.791 

2.50 


.852 

1.211 

1.830 

2.75 


.898 

1.259 

1.872 

3.00 


.932 

1.310 

1.915 

3.25 


.951 

1.355 

1.959 

3.50 



1.400 

2.007 

3.75 



1.441 

2.052 

4.00 



1.482 

2.099 

4.25 



1.512 

2.146 

4.50 



1.540 

2.189 

4.75 



1.551 

2.237 

5.00 



1.540 

2.278 
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TABLE VI. - ELECTRODE COORDINATES FOR RECTILINEAR ION BEAM 
WITH CIRCULAR CROSS SECTION AWAY FROM BEAM EDGE 


z/ a 

r/a 

z/ a 

r/ a 

z / a 

r/a 

z/ a 

r/ a 

X = 

■ 0 

X = 

0.5 

X = 

1.0 

X = 

2.0 

0.2728 

0.962 

0.659 

0.751 

1.000 

0 



.5456 

1.924 

.870 

1.801 

1.173 

1.619 

1.730 

0.994 

.8184 

2.886 

1.108 

2.788 

1.384 

2.667 

1.904 

2.314 

1.0910 

3.848 

1.356 

3.767 

1.611 

3.660 

2.104 

3.400 

1.3640 

4.810 

1.609 

4.733 

1.849 

4.644 

2.317 

4.429 

1.6368 

5.772 

1.873 

5.705 

2.102 

5.625 

2.541 

5.429 

1.9096 

6.734 

2.131 

6 . 666 

2.347 

6.590 

2.777 

6.425 

2.1824 

7.697 

2.393 

7.633 

2.599 

7.564 

3.010 

7.409 

2.4550 

8.659 

2.660 

8.599 

2.857 

8.536 

3.251 

8.391 

2.7280 

9.621 

2.929 

9.561 

3.120 

9.505 

3.504 

9.371 


TABLE VII 


INNER ELECTRODE COORDINATES 


FOR HOLLOW CYLINDRICAL ION BEAM 


z/b 

r/b 


X = 0 

X = 0.5 

X = 1.0 

X = 2.0 

0 

a l. 00000 

b l. 00000 

c l. 00000 

d l. 00000 

.1 

.57984 

.94153 

1.52880 

4.03090 

.2 

.50662 

1.09360 

2.36070 

11.00100 

.3 

.40454 

1.10880 

3.03930 

22.83400 

.4 

.29523 

1.00150 

3.39730 

39.09300 

.5 

.19691 

.81254 

3.35290 

57.09100 

.6 

.12003 

.59487 

2.94820 

72.41500 

.7 

.06687 

.39410 

2.32260 

80.67000 

.8 

.03405 

.23669 

1.64540 

79.51700 

.9 

.01584 

.12904 

1.05100 

69.72300 

1.2 

.00093 

.01181 

.14981 

24.12400 

1.4 

.00009 

.00150 

.02501 

6.98090 

1.6 

0 

.00013 

.00284 

1.33890 

1.8 

0 

.00001 

.00022 

.17162 

0 

.60653 

.60653 

.60653 

.60653 

1.0 

.00674 

.06393 

.60653 

.54598 

2.0 

0 

0 

.00001 

.01480 

3.0 

0 

0 

0 

0 

4.0 




0 


a At beam edge z/b = 0 (see fig. 24(b)). 

^At beam edge z/b = 0.5946 (see 
fig. 24(b)). 

c At beam edge z/b = 1.000 (see 
fig. 24(b)). 

d At beam edge z/b = 1.6818 (see 
fig. 24(b)). 


TABLE VIII. - ELECTRODE COORDINATES FOR CONVERGING CYLINDRICAL FLOW (r/r 0 < l) 


CD 

O 


Y 


X = 


0. 

1 . 732 7307E-02 
3.AA03A09E-02 
5. 1229A39E-02 
6.78065 1 0E-02 
8.A135733E-02 
1 .0021819E-01 
L. 1605A96E-01 
1.316A710E-01 
l . A699567E— 01 
1.6210168E-01 
2.3A02917E-01 
3 . 000A260E-0 1 
3 . 6026066 E-0 1 
A . 1 A79 772 E-01 
A. 6376501 E-01 
5.0727155E-01 

5. A5A2555E-01 
5.7833529E-01 
6.061 1091 E-0 1 
6.288652AE-01 

6. A671565E-01 
6. 5978537E— 01 
6.6820A69E-01 
6.7211 226E— 0 1 
6.716561 IE-01 
6.6699A37E-01 
6.5829513E-01 
6 . A5 7 359 3E-01 
6 . 2950262 E-0 1 
6. 09786A0E-01 
5 . 8678002 E-0 1 
5.606 7 1 87 E-0 1 
5.3163729E-01 
A.9982693E-01 
A . 65 35 1 1 A E-0 1 


X 


0 


1.A901 161E-08 
7.3193A61E-03 
1 .A815383E-02 
2 . 2A8 A l 76E-02 
3.032185 l E-02 
3.832A289E-02 
A.6A8770AE-02 
5.A80830AE-02 
6. 3282333E-02 

7 . 1 90596 3E-02 
8.067578BE-02 
1 . 26591 3AE-01 
1.756398AE-01 
2.27A1862E-01 
2.8 1 5 A937E-01 
3. 376 7 502E-0 1 
3.95A57A1 E-0 l 
A.5A57A56E-01 
5. 1A71918E-01 
5.7559668E-01 
6.3692A59E-01 
6.98A3165E-01 
7. 5985788E-01 
8 . 2095 520E-0 1 
8.81 A 8860E-0 1 
9. A 12 3783E-01 
1 . 000000 IE 00 
1 . 0575935E 00 
1 . I 138608E 00 
1 . 16867A9E 00 
1.2219AA3E 00 
1.2 736202E 00 
1 . 32370A2E 00 
l . 3722567E 00 
l . A 1 9 A0A8E 00 
I.A65351AE 00 


Y 


X = 


0. 

1 . 5676782E-02 
3.1 308669E-02 
A.6855A80E-02 
6.228AAA0E-02 
7. 7570A59E-02 
9.269A856E-02 
1 . 076 A 398E-0 1 
1.22A0783E-01 
1 . 3697898E-0 1 
1.5135 19AE-01 
2.2012883E-01 
2 . 8 36 1 7AAE-0 1 
3.A176071 E-0 l 
3 . 9A5 766 1 E-0 1 
A. A2 1 1 587E-01 
A.8AAA563E-01 
5 . 2 16 A 3A8 E-0 1 
5. 5379397E-0 l 
5.809883AE-01 
6.0332389E-01 
6.2090A67E-01 
6 . 3 38 A 26 3E-0 1 
6 . A22580A E-0 1 
6.A628109E-01 
6 . A6052 1 7E-0 1 
6.A172279E-01 
6 . 33 A5532 E-0 1 
6.2 1A230AE-01 
6 . 05 808A0E-0 1 
5. 8680097E-01 
5.6A593A2E-01 
5.3937621 E-0 1 
5.11 329A2 E-0 l 
A. 8061 167E-01 
A. A73AA51E-01 


X 


0.05 


1.0 1 3 A 75AE-01 
1.0187791E-01 
1.03A3715E-01 
1 . 059AA28E-01 
1.0929112E-01 
1. 1336571E-01 
1. 18067A7E-01 
1.2331023E-01 
1.2902313E-01 
1 . 35 1 A885E-0 1 
l . A 16A060E-01 
1.78A6812E-01 
2.2076759E-01 
2 . 6709 1 89E-0 l 
3. 1657337E-01 
3.6859320E-01 
A. 2266027E-01 
A.7835662E-01 
5.35311 78E-0 1 
5.9318770E-01 
6.5167089E-01 
7. 10A6750E-01 
7.6930026E-01 
8. 2790759E-01 
8 . 860 A 3 32E-0 1 
9.A3A7788E-01 
1 . 000000 1 E 00 
1 . 055A203E 00 
1.1095738E 00 
1 . 1623260E 00 
1.2135783E 00 
1.26327 AAE 00 
1.31 1 A092E 00 
1 . 3580369E 00 
l . AO 32 8 1 3E 00 
1.AA73AA8E 00 


Y 


X = 


0 . 

1 . A55920AE-02 
2.9092681E-02 
A. 35757A2E-02 
5.7985585E-02 
7.2301837E-02 
8.6506767E-02 
1 . 0058527 E—0 1 
1 . 1 A52 A62E-01 
1.2831A23E-01 
1 . A19A529E-01 
2.0750718E-01 
2.68A1175E-01 
3.2AA3286E-01 
3. 75A8927E-01 
A. 2 1 56789E-01 
A.6269212E-01 
A . 9890723E-0 1 
5. 302 7326E-01 
5.5686 167E-0 l 
5.7875361E-01 
5.9603973E-01 
6.0881 998E-0 l 
6 . 1 720A67E— 0 1 
6.213 l A92 E-0 1 
6.212 8 33AE-0 1 
6.1 725A39E— 01 
6 . 0938 A72E-0 l 
5.978A291E-01 
5 . 82 8076A E-0 1 
5 . 6AA6637 E-0 1 
5. A301 158E-01 
5. 186353AE-01 
A.91522A0E-01 
A. 6183891 E-0 1 
A. 2971 812E-01 


X 


0.1 


1.6557 IAAE-0 l 
1 . 6590387E-01 
1.6689A7AE-01 
1.68525A1E-01 
1 . 7076761E-01 
l . 7358633E-0 1 
1.769A319E-01 
1 . 8079878E-0 1 
1 . 851 1509E-01 
l . 8985595E-01 
1 .9A9883AE-01 
2.2557276E-01 
2 . 625AA88E-01 
3.0A25158E-01 
3.A963360E-01 
3.979387AE-01 
A. A858510E-01 
5. 0109286E-01 
5 . 550 A799E-0 l 
6. 1008136E-01 
6.6585655E-01 
7.2206221E-01 
7.78A0775E-01 
8. 3A62061E-01 
8.90AA553E-01 
9 . A56AA87E-01 
1 . 0000001 E 00 
1.05331AAE 00 
1.105A160E 00 
1 . 1561630E 00 
1.205AA91E 00 
1.2532 106E 00 
1.299A3AAE 00 
1.3AA1679E 00 
1 . 3875291E 00 
1.A29718AE 00 


TABLE VIII. - CONTINUED. ELECTRODE COORDINATES FOR CONVERGING CYLINDRICAL FLOW (r/r Q < l) 


Y 

X 

Y 

X 

Y 

X 

X = 

0.2 

X = ( 

0.4 

X = 

0.6 


0. 

1.2823116E-02 
2.5632195E-02 
3.8413405E-02 
5 • 1 1 5 3300E-02 
6 • 383898 3E-02 
7.6458250E-02 
8 • 8999655 E— 02 
1 • 0145261 E— 0 1 
1.1380736E-01 
1 .2605502E-01 
1.8543078E-01 
2.4121209E-01 
2 • 9296916E-01 
3 • 4046437 E-0 1 
3*8357151 E-01 
4.2223080E-01 
4 • 5642468E-01 
4 • 86 1 6425 E-0 1 
5*1 148 127E-01 
5.3242416E-01 
5 • 490558 7 E-01 
5*6145254E-01 
5.6970347E-01 
5 • 739 1 128E-01 
5.7419226E-01 
5 • 7067636E-0 1 
5.6350801E-01 
5.5284489E-01 
5 • 3885 795 E-0 1 
5*2 1 7291 3E-0 1 
5 . 0164882E-01 
4.7881 162E-0 1 
4.5340990E-01 
4*2562461 E-01 
3 *956 1 192E-01 


2 • 65155 14E-0 1 
2.6536427E-01 
2 • 6599005E-0 1 
2.6702865E-01 
2 • 6847383E-0 1 
2.7031718E-01 
2.7 2 54844E-0 1 
2 • 75 1 5606E-0 1 
2.7812724E-01 
2.8144868E-01 
2.85 1 0652E-0 1 
3.0796298E-01 
3.3727531E-01 
3*71 72569E-0 1 
4. 1029850E-01 
4* 5220320E-0 1 
4.9680501E-01 
5 • 4357539E-01 
5.9205984E-01 
6.4185706E-01 
6.9260486E-01 
7.4397113E-01 
7.9564809E-01 
8.4734851 E-0 l 
8 • 9880406E-0 1 
9.4976479E-01 
1 • 0000001 E 00 
1.049300 7 E 00 
1.0974816E 00 
1.1443867E 00 
1.1898940E 00 
1 . 2339228E 00 
1.2764422E 00 
1.3174812E 00 
1.3571417E 00 
1.3956122E 00 


0. 

1.0285877E-02 
2. 0564585 E- 02 
3.0828995E-02 
4. 1072036E-02 
5 • 1286758 E-02 
6. 1466322E-02 
7.1 604072 E-02 
8.1 69 3500 E— 02 
9. 1728344E-02 
1.0 1 70255E— 0 1 
1 • 504689 1 E-0 1 
1 • 9695296 E— 01 
2. 406 7456E-01 
2.8128519E-01 
3. 1 853940E-01 
3.5226841E-01 
3.8236013E-01 
4.0874477E-01 
4.3138507E-01 
4.502696 1 E-01 
4.6540889E-01 
4.7683284E-01 
4. 8458929E-01 
4. 8874348E-01 
4.8937808E-01 
4.8659261E-01 
4 • 805041 5 E-0 1 
4. 7124676E— 01 
4.5897142 E-01 
4. 4384486 E-01 
4.2604782E-01 
4.0577304E-01 
3.8322098E-01 
3 • 5 859436 E-0 1 
3 • 3208930E-01 


4. 1059425E-01 
4 • 1072260E— 0 1 
4.111074 1 E-0 1 
4. I 174780E-01 
4. 1264259E-0 1 
4. I 378975E-0 1 
4. 151 8709E-01 
4. 168317 IE-01 
4. 1 872061E-01 
4. 208501 IE-01 
4.232 1 626E-0 1 
4 • 3844246E-0 1 
4.5887626E-01 
4. 8387 1 80E-01 
5. 1279983E-01 
5.4507865E-01 
5.80 1 8293E— 0 1 
6.1 76 3994E-0 1 
6 • 5702249E-0 1 
6 • 9 794099E-0 1 
7 • 4003679E-01 
7 • 8297635E-0 1 
8 • 26447 1 2E-0 1 
8 • 701 5477E-0 1 
9 • 1 382 142E-0 1 
9.571 8505E-0 1 
1 . 000000 1 E 00 
1 . 0420388E 00 
1 • 0830936E 00 
1 • 122981 1 E 00 
1.1615464E 00 
1 • 1986695E 00 
1.2342733E 00 
1 • 2683344E 00 
1 . 3008966E 00 
1.3320895E 00 


0. 

8.4177403E— 03 
1 . 6830868E— 02 
2. 5234791 E— 02 
3.3624927E-02 
4.1 996744E— 02 
5 • 0345736E— 02 
5.8667475E-02 
6. 69575 79E-02 
7.521 1752E— 02 
8 • 342 5 789E-02 
1 • 2375295E— 01 
1.6250279E-01 
1 • 992 7 1 34E-01 
2.3372754E-01 
2 • 6560768E-0 1 
2 • 947060 1 E-01 
3.2086524E-01 
3 • 4396825E-01 
3 • 639 3096E— 01 
3. 8069745E-01 
3.942 36 1 5E-01 
4. 045 3754E-01 
4. 1161261E-01 
4. 1549208E— 01 
4. 1622600E— 01 
4. 1 388 389E— 0 1 
4. 085 5455E-01 
4.0034621 E— 0 1 
3.8938621E-01 
3 • 7582074E— 01 
3 • 598 1 377E-01 
3. 4154598E-01 
3.2121 273E— 01 
2 • 9902 162 E-01 
2.751 8806E-01 


5.1 765 506E-0 1 
5.1774R04E-01 
5 • 1 802705E-0 1 
5 . 18491 55E-0 1 
5.1914 120E-0 1 
5. 1997507E-01 
5 . 20992 38E-0 1 
5.221918 1 E-0 1 
5.23572 1 3E-0 1 
5.2513174E-01 
5 • 2686889E-0 1 
5.38 14775E-01 
5 • 5352730E-01 
5.7266 1 29E-0 1 
5.95 17208E-0 1 
6. 2067297E-01 
6.4878311 E-0 1 
6.791 3482E-0 1 
7.11 37643E-0 1 
7.4517 286E-0 1 
7.8020428E-01 
8 • 16 16473E-0 1 
8. 5276045E-01 
8. 8970881E-01 
9.2673761E-01 
9.6358502E-01 
1 • 000000 1 E 00 
1 • 0357440E 00 
1.0705919E 00 
1 . 1 043358E 00 
1 . 1 367 877E 00 
1.1677 840E 00 
1.1971918E 00 
1.22491 57E 00 
1 • 2509091 E 00 
1.2751 882E 00 


62 


TABLE VIII. - CONCLUDED. ELECTRODE COORDINATES FOR 
CONVERGING CYLINDRICAL FLOW (r/r Q < l) 


Y 

X 

Y 

X 

X = 

0.8 

X = 

1.0 


0 . 

6. 9609271 E- 03 
1.391 860 3 E— 02 
2.0869779 E— 02 
2.7811222E-02 
3.473971 IE-02 
4. 1652042E-02 
4.8545030E-02 
5 • 541 555 3E— 02 
6. 2260466 E— 02 
6.907672 IE— 02 
1.0262334E-01 
1.3501 944E— 01 
1. 6594096 E— 01 
1 . 95 1 0403E— 0 1 
2 • 2226539E— 01 
2.4722046E— 01 
2 . 6979982E— 01 
2.8986560E— 01 
3.073081 IE-01 
3. 2204290 E— 01 
3. 3400864 E— 01 
3 .431657 1 E— 01 
3. 4949506 E— 01 
3. 5299808 E- 01 
3. 5369632E-01 
3.5163151 E-0 1 
3.4686594E— Ol 
3 . 39482 1 3E—0 1 
3.2958314E— 01 
3. 1729210E-01 
3.0275180E-01 
2.8612359E-01 
2 • 6758647E-0 1 
2 . 4733526E— 01 
2 . 2557906E-01 


6.0113725E-01 
6.0 1208 79E-0 1 
6.0142 340 E-0 1 
6.0178094E-01 
6.0228104E-01 
6.0292333E-01 
6 . 0370734E-0 1 
6 . 046 3250E-0 1 
6.0569791E— 01 
6.0690294E-01 
6.0824649E— 01 
6.1700452E-01 
6 . 290 37 1 8E — 0 1 
6.4413850E— 01 
6.6206996E— 01 
6.825721 7E-0 1 
7.0537417E-01 
7.3020023E-01 
7.5677392E— 01 
7.8482060E-01 
8. 1406850E— 01 
8.4424926E— 01 
8. 750979 4E— 01 
9.0635314E— 01 
9.3775697E— 01 
9.6905564E-01 
1 . 000000 1 E 00 
1.0303469E 00 
1 . 059 8599E 00 
1.0883116E 00 
1.1154850E 00 
1. 1411756E 00 
1 . 1651936E 00 
1.1873662E 00 
1 . 2075390E 00 
1.2255792E 00 


0 . 

5. 7926631 E — 03 
1. 1582931E-02 
1. 7368407 E— 02 
2.3146702E-02 
2. 8915439 E— 02 
3.4672235E-02 
4.0414739E-02 
4.6140609E— 02 
5. 1847492E-02 
5.7533106E-02 
8.5562512E-02 
l . 1272692E— 01 
1 . 3876833E— 01 
1.6345069E-01 
l . 8656168E-01 
2.0791291 E-0 1 
2.2733941E-01 
2.4469884E— 01 
2 . 5987023E-01 
2.7275310E-01 
2.8326694E-01 
2.9135070E— 01 
2 . 969627 3 E— 01 
3.000812 IE— 01 
3.0070402E-01 
2 .98 84966 E-0 1 
2.9455707E-01 
2.8788637E— 01 
2,7891851 E— 0 1 
2.6775547E-01 
2.5451967E-01 
2. 3935304E-01 
2. 2 24 l 644 E— 01 
2.0388784E-01 
1 . 8396113E-01 


6.6808197E— 01 
6.681 3859E— Ol 
6.6830828E-01 
6 . 6859 l 09E— 01 
6.6898677E— Ol 
6.6949506E-01 
6.7011 5 76E-0 1 
6. 7084841E— 01 
6.7 169255E-01 
6.7264785E-01 
6.7371358E-01 
6.8067639E-01 
6.9028537E-01 
7.0241039E-01 
7. 1689513E-01 
7 . 3 3562 34 E— 01 
7.522 1908E— 01 
7.7266063E-01 
7 . 9467 334E-0 1 
8.180 3695E-01 
8. 4252597E-01 
8 . 679 1050E-0 1 
8 . 9395707E-01 
9. 204291 2E-0 1 
9. 4708763E— 01 
9. 7369186E-01 
1 . 000000 l E 00 
l • 0257705E 00 
1 . 0507622E 00 
1 . 074736 1 E 00 
l . 0974551 E 00 
1.1 1 86845E 00 
1.1381 906E 00 
l . 1557378E 00 
1 . 1710823E 00 
1 . 1839612E 00 


TABLE EX. - ELECTRODE COORDINATES FOR DIVERGING CYLINDRICAL FLOW (r/r 0 > 1) 


Y 

X 

Y 

X 

Y 

X 

X 

= 0 

X = 

0.2 

X = 

0.4 


0 . 

1. 7579666E-02 
3.5412857 E-02 
5 . 3500735 E-02 
7. 1 844539E-02 
9.0445459E-02 
1 . 0930475E-01 
1.2842 366E-0 1 
1 . 4780346E-01 
1.6744551 E-0 1 
1.87351 06E-0 1 
2 . 90880 1 OE-O 1 
4.0121 142E-01 
5. 1854261 E-0 1 
6 . 4309 365 E-0 1 
7. 751 1073E-01 
9. 1486998 E-0 1 
l . 062681 3E 00 
1.2 188907E 00 
1.3838814E 00 
1.5580718E 00 
1.741 9094E 00 
1.9358589E 00 
2 . 1403832E 00 
2.3559130E 00 
2 . 582 8076E 00 
2 . 82 l 3036E 00 
3.0714571E 00 
3 • 3330855E 00 
3 . 605 7 34 7 E 00 
3 . 8887087E 00 
4.18 1 2666E 00 
4.4832159E 00 
4.7965171E 00 
5.1 298906E 00 
5.5 149460E 00 


-7.4505806E-09 
7 . 1 382 374E-03 
1.4091402E-02 
2 . 0854697E-02 
2. 7424648 E-02 
3 . 3796206E-02 
3.9965063E-02 
4.5926630E-02 
5. 1676 139E-02 
5.7209 1 79E— 02 
6.2520504E-02 
8.5579171E-02 

1 . 02 3 1 9 1 5E— 0 1 
1. 1201806E-01 
1. 1387070E-01 
1 . 0697266E-0 1 
9 . 0299 398E-02 
6.2681 124E-02 

2.2 770494E-02 
-3.0996077E-02 
-1 .0044712E-01 

- 1 . 877 3461 E-0 1 
-2.9540533E-01 
-4.2648645E-01 
-5.8458941E-01 
-7. 7403412E-01 
- 1 . 0000006E 00 
-1.2687183E 00 

— 1.587713 7E 00 
-1.9661546E 00 
-2.4153753E 00 
-2.9497579E 00 
-3.5883874E 00 
-4. 3585593E 00 
-5.3044917E 00 
-6. 5 149492E 00 


0 . 

2.3334026E-02 
4.6690750E-02 
7.0092498E— 02 
9 .356091 1 E-02 
1.171 1662E-01 
1.4077909E-01 
1 .6456644E-01 
1 .8849541E-01 
2. 1258137E-01 
2 . 3683830E-0 1 
3 . 6 109524E-0 l 
4.91 1 22 76 E-0 1 
6 . 277 3326 E-0 1 
7.7 152793E-01 
9.2301402E-01 
1 .082665 1 E 00 
1.2509553E 00 
1.4283786E 00 
1.61 5461 6E 00 
1.8127676E 00 
2.0209011E 00 
2 • 240501 2 E 00 
2.4722355E 00 
2.7167856E 00 
2 • 9748363E 00 
3.2470645E 00 
3.5341515E 00 
3.8368448E 00 
4. 156 1400E 00 
4.4937255E 00 
4. 8530455E 00 
5.2419587E 00 
5.6802182E 00 
6 . 2265034E 00 


3 . 3672 364E-0 1 
3 . 3680530E-0 1 

3.3 704 796E-0 1 
3. 3744437E-01 
3. 3798319E-01 
3 . 3864893E— 0 1 
3 . 3942334E-0 1 
3. 4028590E-01 
3.412 1434E-0 1 
3.42 1 8574E-01 
3.4317659E-01 
3.4762563E-01 
3.4934852E-01 
3.46178 l 4E-0 1 
3. 36 32539E-0 1 
3. 1820042E-01 
2 . 902698 l E-0 l 
2 . 5095502E-01 
1. 9855171 E-01 
1.31 15807E-01 
4.6601564E-02 
-5.7638660E-02 
-1.8452465E-01 
-3.3756696E-01 
-5.2095791E-01 
-7. 397361 3E-01 
- 1 . 0000007E 00 
-1.3091990E 00 
-1.6765402E 00 
-2.1136354E 00 
-2 . 6355836E 00 
-3.2630137E 00 
-4. 026448 IE 00 
-4.9773338E 00 
-6.2246597E 00 


0 . 

2.8 1 56436E-02 
5.6329107 E-02 
8.45341 38E-02 
1. 1 278749E-01 
1 .41 10482E-01 
1.6950149E-01 
1 . 9799242E-01 
2 . 2659205E-01 
2 . 553 1434E-01 
2.841 7275E-01 
4. 3093650E-01 
5.8280022E-01 
7 .408705 1 E-01 
9. 0604094 E— 01 
1 . 0790749E OO 
l • 2606677E 00 
1.4514924E 00 
1.6522315E 00 
1 . 8635987E 00 
2.0863557E 00 
2.321 3259E 00 
2. 5694041 E 00 
2.8315707E 00 
3. 1089087E 00 
3 . 4026390E 00 
3.7 14 1906E 00 
4.0453506E 00 
4.3985717E 00 
4. 7776297E 00 
5.1890756E 00 
5.6457168E 00 
6. 1 762 1 1 3 E 00 
6.8593919E 00 


6. 1 308981E-01 
6. 1 308096E-01 
6. 1305435E-01 
6.1 300728E-0 1 
6 . 1 293600E-0 1 
6. 1283527E-0 1 
6. 1269876E-01 
6. 1251867E-01 
6. 1228603E-01 
6. 1199106E-01 
6.1 1 62 354E-0 1 
6 . 082 7674E-0 1 
6.0123022E-01 
5.8880171E-01 
5.6930873E-01 
5.4107837E-01 
5.0240499E-01 
4.5 149205E-0 1 
3 . 8638656E-0 l 
3 . 0490550E-0 1 
2 . 0455 766E-0 1 
8.245 1 626E-02 
-6.4813845E-02 
-2 .4128 343E-0 1 
-4.5181609E-01 
-7.0230833E-01 
-1.0000008E 00 
-1.3539232C 00 
— 1.7755 879E 00 
-2 .280 1633E 00 
-2 . 8886706E 00 
-3.6326426E 00 
-4.5658393E 00 
-5 . 80 30006E 00 


CD 

W 


TABLE IX. - CONTINUED . ELECTRODE COORDINATES FOR DIVERGING CYLINDRICAL FLOW (r/r 0 > l) 


CD 

v^ 


Y 

X 

Y 

X 

Y 

X 

X = 

0.6 

X = O.i 

S7825 

X = 

0.8 


0. 

3 • 2985 843E-02 
6.5985871E-02 
9 .9014231E-02 
1.3208498E-01 
1.652 1 205 E-0 1 
1 • 984092 1 E— 0 1 
2.3 1 69004 E-01 
2 • 6506783E-01 
2.9855577E-01 
3.3216655E-01 
5 • 0249035 E-01 
6. 7762592E-01 
8.58821 40E-01 
1 . 047 1626E 00 
1 • 2436208E 00 
1 • 449 1020E 00 
1 • 6644896 E 00 
1 • 89068 1 OE 00 
2.1286179E 00 
2.3793153E 00 
2 • 6438909E 00 
2 • 9235959E 00 
3.2 1 98635E 00 
3. 534378 l E 00 
3 • 869 1 927 E 00 
4.22 69456E 00 
4.61 12645E 00 
5 • 0275726E 00 
5. 4847822E 00 
5 . 9992 197 E 00 
6 • 6050682E 00 
7 • 390061 3E 00 
8 • 7077599E 00 


8 • 898 1 370E-0 1 
8 • 8975 743E-0 1 
8 • 895 8 788E-0 1 
8.8930384E-01 
8 • 8890299E-0 1 
8.883821 7E-0 1 
8.8773733E-01 
8 • 8696 369E-0 1 
8.8605536E-01 
8.8500668E-01 
8.8380989E-01 
8.753 1 930E-0 l 
8.6176 1 67E-0 1 
8 • 41 748 1 8E-0 1 
8.1 37 3854E-0 1 
7.760742 3E-0 1 
7.2697217E-01 
6.6448924E-01 
5.8646935E-01 
4. 904 7 390E-0 1 
3 • 7 369 794E-0 1 
2.3286609E-01 
6.4101368E-02 
-1.3724078E-01 
-3.7679445E-01 
-6.6149026E-01 
- 1 . 0000009E 00 
-1.4034344E 00 
- 1 . 8864979E 00 
-2.4696444E 00 
-3. 1835391E 00 
—4. 0799958E 00 
-5.2666563E 00 
-7. 0972426E 00 


0. 

3.4908992E-02 
6.983 1 797E-02 
1.0478216E-01 
1.3977382E-01 
1.7482034E-01 
2 • 0993522E-01 
2.4513177E-01 
2 • 8042 31 OE-0 1 
3. 1582220E-01 
3.51 34 164E-01 
5.311 7205 E-0 1 
7. 1576786E-01 
9.0642015E-01 
1.1 04277 1 E 00 
1 • 3103836E 00 
l • 5257256E 00 
1.7512717E 00 
1. 9880102 E 00 
2.2 369844E 00 
2.4993252E 00 
2. 7762878E 00 
3.0692970E 0.0 
3. 3800088E 00 
3. 7 104 1 1 7E 00 
4. 062989 7 E 00 
4 • 44 1 0203E 00 
4 • 849 1 306 E 00 
5.2944126E 00 
5 • 788802 3E 00 
6 • 354770 1 E 00 
7 • 04 1 4249E 00 
7 . 9866989E 00 


1 • 000005 8E 00 
9. 999 3455E-0 1 
9.9972075E-01 
9.9936248E-01 
9.9885848E-01 
9.9820534E-01 
9 • 9739990E-0 1 
9.9643780E-01 
9.9531 379E-0 1 
9 . 940226 1 E-0 1 
9.9255724E-01 
9 • 8236085E-0 1 
9.6655575E-01 
9.4382402E-01 
9. 1266372E-01 
8.7142142E-01 
8. 1828870E-01 
7.5127 129E-0 1 
6.681 3828E-0 1 
5 • 663529 3E-0 1 
4.4298564E-01 
2 • 94603 78E-0 1 
1. 1713217E-01 
-9.4329997E-02 
-3.4575500E-01 
-6.4453526E-01 
— 1 . 0000009E 00 
- 1 . 4242450E 00 
- 1 • 9335490E 00 
-2.5511064E 00 
-3.3129490E 00 
-4.2834724E 00 
-5 • 6 1 1 1 254E 00 


0. 

3.7951 939E-02 
7 • 59 1 7358 E-02 
1. 1390967E-01 
1 • 519423 1 E-01 
1 .9002859E— 01 
2.2818 174E-01 
2 • 664 1 490E-0 1 
3.0474097E-01 
3.431 7285E-0 1 
3.81723 1 7E-0 1 
5.7668148E-01 
7.7638677E-01 
9 • 82 1 9366E-0 1 
1.1953417E 00 
1 • 4169867E 00 
1 • 6482389E 00 
1.8902038E 00 
2 • 144021 2E 00 
2 . 4109060E 00 
2.6921904E 00 
2.9893748E 00 
3.3041949E 00 
3.6387179E 00 
3 • 9954974E 00 
4. 3778323E 00 
4 • 7902 31 9E 00 
5.2392960E 00 
5 • 735501 8E 00 
6.2971761E 00 
6.9606043E 00 
7.8122556E 00 


1 . 1 743571E 00 
1 • 1742649E 00 
1 • 1 739884E 00 
1 • 1735258E 00 
1.1728760E 00 
l . 1 720365E 00 
1.17 10039E 00 
1.1697749E 00 
1 • 1 68 3443E 00 
1. 166 7078E 00 
1.1648594E 00 
1 • 1522043E 00 
1.133 1048E 00 
l • 1063208E 00 
1 • 0703923E 00 
1 • 02 36664E 00 
9 • 642942 3E-0 1 
8 • 9020342E-0 1 
7 • 9904696E— 0 1 
6.881 34 10E-0 1 
5 • 5433637E-01 
3. 9396 785E-0 1 
2 • 0262799E-0 1 
-2.5009155E-02 
-2.9548675E-01 
-6. 1699C17E-01 
- 1 • 00000 10E 00 
-1.4583801E 00 
-2.0113251E 00 
-2.6873249E 00 
-3. 5334548E 00 
-4.6429169E 00 


TABLE DC. - CONCLUDED. ELECTRODE COORDINATES FOR DIVERGING CYLINDRICAL FLOW (r/r Q > l) 


Y 

X 

Y 

X 

Y 

X 

X = 

1.0 


x = : 

L .5 


X = 

2.0 


0 . 

1.4695159E 

00 

0 . 

2. 2609916E 

00 

0 . 

3. 1333990E 

00 

A. 3103422E-02 

1.4693931E 

00 

5.6917496E-02 

2 • 2608008E 

00 

7.2144340E-02 

3. 1331460E 

00 

8 . 6220226E— 02 

1.4690239E 

00 

1. 1384969E-01 

2 • 2602280E 

00 

1 .4430612E-01 

3. 1323860E 

00 

1 . 2936380E— 01 

1 • 4684080E 

00 

1.7081 126E-0 1 

2.2592726E 

00 

2.1 650280E-01 

3.1311191E 

00 

L.7254747E-01 

1 . 4675435E 

00 

2.2781 690E-01 

2.2579330E 

00 

2.8875 186E-0 1 

3 • 1293434E 

00 

2. 1578455E-01 

1 • 4664284E 

00 

2.8488127E-01 

2 • 2562073E 

00 

3.6107073E-01 

3. 1270568E 

00 

2. 5908829E— 01 

1 • 465060 1 E 

00 

3. 4201905E— 01 

2 • 2540934E 

00 

4. 3347694E-0 1 

3. 1242570E 

00 

3.0247179E-01 

1.4634347E 

00 

3. 9924485E— 01 

2.2515879E 

00 

5.0598798E-01 

3 • 1209409E 

00 

3 • 459482 1 E-01 

1.4615491E 

00 

4. 5657336E— 01 

2.2486878E 

00 

5.7862 142E-01 

3. 1171048E 

00 

3 . 895305 1 E— 0 1 

1.4593985E 

00 

5. 1401907E-01 

2.2453883E 

00 

6.5 1 39486E— 01 

3. 1 127447E 

00 

4.3323155E-0I 

1.4569779E 

00 

5. 7159666E— 01 

2 • 24 1 6853E 

00 

7 • 2432590E-0 1 

3. 1078557E 

00 

6.5 396060E-01 

1.4406139E 

00 

8.6197 140E— 01 

2.2169207E 

00 

1.0919652E 00 

3 • 0752689E 

00 

8. 7951779E-01 

1.4164550E 

00 

1 . 1 578182E 00 

2. 1810791E 

00 

1.4662235E 00 

3.02841 54E 

00 

1.11 13659E 00 

1 • 383331 5E 

00 

1 • 4609435 E 00 

2 • 1330031E 

00 

1 • 8493993E 00 

2 • 966049 1 E 

00 

1.3508922E 00 

1.3398136E 

00 

1 • 7731 655 E 00 

2.0712122E 

00 

2 • 243 8950E 00 

2 • 8865396E 

00 

1.5994346E 00 

1 • 2842290E 

00 

2 • 0963447E 00 

1 • 9938900E 

00 

2 • 6522758E 00 

2 . 7878418E 

00 

1.8583150E 00 

1.2 146531E 

00 

2 • 4324308E 00 

1.8988577E 

00 

3 • 0773492E 00 

2 • 667441 3E 

00 

2. 1288799E 00 

1 • 1288795E 

00 

2 • 78352 1 9E 00 

1 • 78352 1 3E 

00 

3 • 5222776E 00 

2 • 5222769E 

00 

2.4125475E 00 

1 . 0243673E 

00 

3.151 9463E 00 

1 • 6447964E 

00 

3.9907363E 00 

2 • 3486246E 

00 

2.7 1 08584E 00 

8.9816294E- 

01 

3. 5403725E 00 

1 • 4789948E 

00 

4.487141 1 E 00 

2. 1419291E 

00 

3.0255391E 00 

7.4679531E- 

01 

3 • 95 1 962 7E 00 

1 • 28 1 6660E 

00 

5.0 16992 8E 00 

1 • 896561 3E 

00 

3. 3585833E 00 

5.661 3253E- 

01 

4.3906007E 00 

1 • 0473700E 

00 

5 . 5874 1 3 3E 00 

1.6054527E 

00 

3.712 3697E 00 

3.51191 39E- 

01 

4.8612435E 00 

7.6934704E 

-01 

6.2080304E 00 

1.2595371E 

00 

4. 0898345E 00 

9 • 5867082E- 

02 

5 • 3704993E 00 

4. 3901993E 

-01 

6. 8925108E 00 

8.4684141E 

-01 

4.4947588E 00 

-2.0745360E- 

01 

5 • 9276203E 00 

4.5198262E 

-02 

7.6614143E 00 

3.5091260E 

-01 

4.9322571E 00 

-5 . 6848438E- 

•01 

6 • 5463442 E 00 

-4.2727039E 

-01 

8 • 5479789E 00 

-2.52 1 5045E 

-01 

5.4096707E 00 

- 1 . 000001 IE 

00 

7 • 2485798 E 00 

-1.0000015E 

00 

9.611 3002 E 00 

— 1 • 0000020E 

00 

5.93831 36E 00 

-1.5195364E 

00 

8 . 072631 2E 00 

-1.7062655E 

00 

1 • 097 1956E 01 

- 1 • 9599242E 

00 

6. 5371969E 00 

-2 • 1 526857E 

00 

9. 0947183E 00 

-2.6036464E 

00 




7.2420303E 00 

-2 • 940498 IE 

00 

1.0501963E 01 

-3.81 39954E 

00 




8.1 316562E 00 

-3. 9596832E 

00 







9. 4393312E 00 

-5 • 40 16352E 

00 








TABLE X. - OUTER ELECTRODE COORDINATES FOR CIRCULAR FLOW ION BEAM 


e 

R 

e 

R 

0 

R 

e 

R 

e 

R 

e 

R 

Q 

R 

G 

R 


X = 

0 


X = o. 

16494 


X = 0.: 

39685 

X " 0. 

62996 

X = o. 

82549 

X - 0- 

95483 

o 

1.4286 

29.794 

2.7774 

10.000 

1.4286 

35.783 

2.3843 

20.000 

1.4286 

30.000 

1.4286 

40.000 

1.4286 

50.000 

1.4286 

.57296 

1.4629 

30.940 

2.8126 

10.573 

1.4758 

36.929 

2.4089 

20.573 

1.4527 

30.573 

1.4426 

41.146 

1.4445 

51.146 

1.4358 

1.14590 

1.4967 

32.086 

2.8464 

11.719 

1.5508 

38.075 

2.4323 

21.719 

1.4978 

31.146 

1.4563 

42.292 

1.4594 

52.292 

1.4421 

1.71890 

1.5302 

33.232 

2.8790 

12.865 

1.6141 

39.221 

2.4546 

22.865 

1.5395 

32.292 

1.4825 

44.011 

1.4799 

53.438 

1.4475 

2.29180 

1.5633 

34.377 

2.9103 

14.011 

1.6712 

40.367 

2.4756 

24.011 

1.5785 

33.438 

1.5072 

45.730 

1.4981 

55.157 

1.4540 

3.43770 

1.6284 

35.523 

2.9403 

15.157 

1.7243 

41.513 

2.4955 

25.157 

1.6153 

34.584 

1.5306 

47.448 

1.5142 

56.875 

1.4585 

4.58370 

1.6921 

36.669 

2.9689 

16.303 

1.7744 

42.659 

2.5142 

26.303 

1.6500 

35.730 

1.5527 

49.167 

1.5281 

58.594 

1.4611 

5.72960 

1.7544 

37.815 

2.9963 

17.448 

1.8221 

43.805 

2.5318 

27.448 

1.6830 

36.875 

1.5735 

50.886 

1.5400 

59.740 

1.4618 

6.87550 

1.8154 

38.961 

3.0223 

18.594 

1.8678 

44.950 

2.5481 

28.594 

1.7144 

38.021 

1.5932 

52.605 

1.5497 



8.02140 

1.8751 

40.107 

3.0470 

19.740 

1.9118 

46.096 

2.5632 

29.740 

1.7443 

39.167 

1.6118 

54.324 

1.5574 



9.16730 

1.9335 

41.253 

3.0703 

20.886 

1.9541 

47.242 

2.5772 

30.886 

1.7728 

40.313 

1.6293 

56.043 

1.5631 



10.31300 

1.9907 

42.399 

3.0923 

22.032 

1.9949 

48.388 

2.5900 

32.032 

1.7999 

41.459 

1.6456 

57.762 

1.5667 



11.45900 

2.0467 

43.545 

3.1129 

23.178 

2.0343 

49.534 

2.6015 

33.178 

1.8258 

42.605 

1.6610 

59.481 

1.5684 



12.60500 

2.1014 

44.691 

3.1322 

24.324 

2.0723 

50.680 

2.6119 

34.324 

1.8505 

43.751 

1.6752 





13.75100 

2.1550 

45.837 

3.1501 

25.470 

2.1091 

51.826 

2.6211 

35.470 

1.8740 

44.897 

1.6885 





14.89700 

2.2073 

46.983 

3.1666 

26.616 

2.1445 

52.972 

2.6290 

36.616 

1.8963 

46.043 

1.7008 





16.04300 

2.2585 

48.128 

3.1818 

27.762 

2.1787 

54.118 

2.6358 

37.762 

1.9175 

47.189 

1.7120 





17.18900 

2.3084 

49.274 

3.1956 

28 . 908 

2.2117 

55.264 

2.6414 

38.908 

1.9376 

48.335 

1.7223 





18.33500 

2.3572 

50.420 

3.2080 

30.054 

2.2434 

56.410 

2.6458 

40.054 

1.9565 

49.481 

1.7316 





19.48100 

2.4047 

51.566 

3.2190 

31.199 

2.2740 

57.555 

2.6489 

41.199 

1.9744 

50.626 

1.7399 





20.62600 

2.4510 

52.712 

3.2286 

32.345 

2.3033 

58.701 

2.6509 

42.345 

1.9912 

51.772 

1.7473 





21.77200 

2.4962 

53.858 

3.2368 

33.491 

2.3315 

59.847 

2.6517 

43.491 

2.0068 

52.918 

1.7537 





22.91800 

2.5400 

55.004 

3.2436 

34.637 

2.3585 



44.637 

2.0215 

54.064 

1.7591 





24.06400 

2.5827 

56.150 

3.2490 





45.783 

2.0350 

55.210 

1.7636 





25.21000 

2.6242 

57.296 

3.2530 





46.929 

2.0475 

56.356 

1.7671 





26.35600 

2.6643 

58.442 

3.2556 





48.075 

2.0590 

57.502 

1.7696 





27.50200 

2.7033 

59.588 

3.2569 





49.221 

2.0693 

58.648 

1.7712 





28.64800 

2.7410 







50.367 

2.0787 

59.794 

1.7719 













51.513 

2.0869 















52.659 

2.0942 















53.805 

2.1003 















54.950 

2.1055 















56.096 

2.1096 















57.242 

2.1126 















58.388 

2.1146 















59.534 

2.1155 








TABLE XI 


INNER ELECTRODE COORDINATES FOR CIRCULAR FLOW ION BEAN 


R 



e 


X 

II 

O 

X = 0.16494 

X « 0.39685 

X = 0.62996 

X = 0.82549 

x = 0.95483 

X = 1.0 

X = 1.2856 

X = 2.0408 

0 

15 . 00000 

15.0000 

15.000 

15.000 

15.000 

15.000 

15.000 

15.000 

15.000 

.02 

14.99000 

14.9910 

14.993 

14.994 

14.996 

14.997 

14.997 

14.999 

15.004 

. 05 

14.93800 

14.9460 

14.956 

14.966 

14.975 

14.981 

14.983 

14.995 

15.028 

.10 

14.76000 

14.7900 

14.831 

14.872 

14.906 

14.929 

14.937 

14.988 

15.120 

. 15 

14.47600 

14.5410 

14.634 

14.727 

14.805 

14.856 

14.874 

14.988 

15.285 

.20 

14.09400 

14.2110 

14.375 

14.541 

14.679 

14.771 

14.803 

15 . 005 

15.534 

.25 

13.62400 

13.8070 

14.065 

14.324 

14.541 

14.684 

14.734 

15.051 

15.878 

. 30 

13.07400 

13.3390 

13.711 

14.085 

14.398 

14.606 

14.678 

15.136 

16.330 

. 35 

12.45200 

12.8150 

13.324 

13.835 

14.263 

14.546 

14.645 

15.270 

16 . 903 

.40 

11.76500 

12.2420 

12.911 

13.582 

14.145 

14.517 

14.646 

15.467 

17.612 

.45 

11.01800 

11.6280 

12.482 

13.339 

14.056 

14.529 

14.695 

15.740 

18 . 475 

.50 

10.21700 

10.9800 

12.047 

13.115 

14.008 

14.597 

14.803 

16.104 

19.513 

.52 

9.88320 

10.7130 

11.874 

13.034 

14.003 

14.643 

14.867 

16.279 

19.983 

. 55 

9.36760 

10.3060 

11.616 

12.923 

14.015 

14.736 

14.987 

16.577 

20.752 

.57 

9.01480 

10.0310 

11.447 

12.859 

14.037 

14.815 

15.086 

16.801 

21.311 

.60 

8.47260 

9.6135 

11.200 

12.779 

14.094 

14.962 

15.265 

17.179 

22.224 

.62 

8.10270 

9.3331 

11.041 

12.738 

14.151 

15.082 

15.407 

17.462 

22.886 

. 65 

7.53590 

8.9108 

10.813 

12.698 

14.266 

15.298 

15.658 

17.935 

23.967 

.67 

7 . 15020 

8.6292 

10.670 

12.689 

14.365 

15.469 

15.853 

18.287 

24.750 

.70 

6.56040 

8.2087 

10.473 

12.704 

14.552 

15.768 

16.192 

18.872 

26.027 

.72 

6.15990 

7.9311 

10.355 

12.736 

14.705 

16.000 

16.451 

19.304 

26.954 

. 75 

5 .54860 

7.5213 

10.203 

12.822 

14.982 

16.400 

16.894 

20.021 

28.466 

. 77 

5.13420 

7.2548 

10.121 

12.908 

15.201 

16.706 

17.230 

20.547 

29.564 

.80 

4.50240 

6.8692 

10.034 

13.087 

15.588 

17.226 

17.797 

21.413 

31.364 

.82 

4.07470 

6.6251 

10.004 

13.242 

15.887 

17.618 

18.221 

22.044 

32.678 

.85 

3.42350 

6.2858 

10.009 

13.536 

16.402 

18.277 

18.930 

23.078 

34.851 

.87 

2. 98300 

6.0837 

10.052 

13.776 

16.793 

18.766 

19.453 

23.826 

36.459 

.90 

2.31290 

5.8311 

10.185 

14.208 

17.454 

19.576 

20.316 

25.043 

39.176 

.92 

1 . 86010 

5.7081 

10.325 

14.546 

17.945 

20.169 

20.945 

25.920 

41.254 

. 95 

1.17160 

5.6178 

10.621 

15.131 

18.758 

21.138 

21.970 

27.338 

44.981 

. 97 

. 70657 

5.6381 

10.878 

15.572 

19.352 

21.837 

22.709 

28.356 

48.159 

1.00 

0 

5.8124 

11.352 

16.310 

20.318 

22.968 

23.900 

30.000 

58.000 
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Figure 1. - Space-charge-limited flow between infinite plane electrodes, j 
(See appendix C for rescaling.) 
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Figure 2. - Nomographic solution of i mv^ = q((p Q - cp). (See appendix C for 
rescaling. ) 
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Figure 4. - Converging flow between coaxial cylinders, j/l = 3. 43xlO“ 7 A” 1 / 2 (cp 0 - qp) 3 / 2 ( rp2) _1 f 0 r r < r 0 ; 
note that J = (j/l)/27rr. (See appendix C for rescaling.) 
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Figure 5. - Diverging flow between coaxial cylinders. J/l = 3. 43xlO“ 7 A“ 1 / 2 (q> 0 
j = (j/i)/27rr. (See appendix C for rescaling.) 
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Figure 7. - Converging flow between concentric spheres. J =. 6.87xl0 _7 A" 1 / 8 (<p o - T )5/2 a -2 f or r < r . note that 
J =* J/47rr2, (See appendix C for rescaling.) 
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Figure 8. - Diverging flow between concentric spheres. J = 6 . 87xlO" 7 A _1 / 2 (cp 0 - cp)^/ 2 a“ 2 
for r > r 0 ; note that j = j/4irr 2 . (See appendix C for rescaling.) 
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Figure 9. - Circular flow. 
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Hyperbolic flow. J * 4. 92xlO _7 A _1 / 2 ((p 0 
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Figure 12. - Angle between trajectories and normals to equipotentials for in- 
clined planes (eqs. (42) and (49)). 
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Figure 13. - Space-charge-limited flow between inclined planes 
J = 5.467Xl0- 8 A- 1 /2( (po _ ,p)3/2 r -2j ’jL p(e)J. 

(See appendix C for rescaling.) 
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Figure 14. - Potential distribution between inclined plane 
electrodes (space-charge limited). 
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Figure 15. - General trajectory plot for flow between co- 
axial right circular cones (eq. (59)). 
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Figure 16. - Angle between trajectories and normals to equipotentials 
coaxial right circular conical electrodes (eq. (60)). 
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Figure 17. - Space-charge-limited flow between coaxial cones. 
J 0 = S.^xlO'^- 1 / 2 ^ - <p) 3 / 2 r- 2 S(a 2 - c^). 

(See appendix C for rescaling.) 
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Figure 18. - Space- charge -limited, potential distribution between coaxial 
cones . 
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Figure 19. - Current density as a function of electric field 
at emitter (ref. 15). 
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Figure 20. - Normalized potential distribution for partial 
space-charge flow (ref. 15). 


> 


Normalized electric field, dv/dx 


88 



0 .2 .4 .6 .8 1.0 


Normalized distance from emitter, x 

Figure 21. - Electric field in space between emitter and col- 
lector for several current densities (ref. 15). 





Figure 22. - Typical electrode shapes for a 
rectilinear slab ion beam. 
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Figure 23. - Dominant terms in series expansion of potential 
function near edge of a rectilinear ion beam -with circular 
cross section (eq. (79)). 
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(a) Outer electrode shapes. 

Figure 24. - Typical electrode 

shapes for a rectilinear ion beam 
with circular cross section. 
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(b) Inner electrode shapes. 

Figure 24. - Concluded. Typical electrode shapes for a rectilinear ion beam with circu- 
lar cross section. 
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Figure 26. - Electrode shapes for diverging flow between coaxial cylinders (r/r 0 > l). 


J 


a b 


Figure 27. - Electrode shapes for circular flow. X^ 
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Figure 28. - Electrode shapes for hyperbolic flow. 
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Figure 29. - Approximate ratio of ion beam spreading due to a surface irregularity to ther- 
mal spreading. 
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Figure 30. - Effect of electrode tilt on current den- 
sity in a plane diode (ref. 10). (Rectangular emit- 
ter tilted about axis parallel to side and disk 
emitter. ) 
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Figure 31. - Effect of electrode tilt on current density in a 
plane diode with disk emitter of optimum aspect ratio 
D/S = 1.42i (ref. 10). 
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